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RESEARCH MEMORANDUM 


DAMPING IN PITCH OF LOW-ASPECT-RATIO WINGS 
AT SUBSONIC AND SUPERSONIC SPEEDS 


By Murray Tobak 
SUMMARY 


The concept of indicial functions is applied to the analysis of the 
aerodynamic phenomena associated with the short-period pitching mode of 
wings in subsonic and supersonic flight. Simple physical relationships 
are pointed out and are used to study the effect on the rotary-damping- 
moment coefficient of changes in center-of-gravity position, Mach number, 
aspect ratio, plan form, frequency, and thickness. Qualitative conclu- 
sions are drawn from the results of this investigation and are compared 
with the results of experiments for a series of low-aspect-ratio wing- 
body combinations having triangular, swept, and unswept wing plan forms. 


Results of the experimental investigation, which were obtained by a 
single-degree-of-freedom free-oscillation technique over the Mach number 
ranges 0.6 to 0.9 and 1.2 to 1.9, were in good agreement with the results 
of theoretical computations. The predictions of ranges of supersonic 
Mach numLer and center-of-gravity positions over which dynamic instabil- 
ity may be expected, of the beneficial effect on the damping moment of a 
reduction in aspect ratio, and of only a small effect of thickness on the 
damping moment were confirmed by the experimental results. 


The occurrence at high subsonic Mach numbers of small-amplitude 
self-sustained pitching oscillations is noted, and a hypothesis is 
advanced for its explanation. 


INTRODUCTION 


In the classical study of the longitudinal motion of an aircraft, 
it is usually found that the motion resulting from a small equilibrium- 
destroying disturbance consists of two modes: one, & lightly damped, low- 
frequency motion at essentially constant attitude, called the phugoid 
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oscillation; the other, a rotary-pitching and plunging oscillation of 
high frequency (relative to the phugoid frequency) called the short- 
period oscillation. The phugoid oscillation has generally been described 
as resulting from a slow interchange of potential and kinetic energy as 
the aircraft experiences periodic variations in airspeed and altitude. 
The character of the phugoid motion as influenced by airspeed, altitude, 
and aircraft geometry has been well understood for some time (see, e.g., 
ref. 1). The short-period motion, on the other hand, having in the past 
been found to be highly damped and of short duration, has been the cause 
of no concern. Its characteristics therefore have not been as fully 
investigated as those of the phugoid oscillation. With the advent of 
flight at speeds approaching and exceeding the speed of sound, however, 
the loss of rotary damping occurring in practically all aircraft at 
speeds near the sonic speed has caused renewed interest in the short- 
period pitching mode. Unlike the easily controlled phugoid oscillation, 
the deterioration of damping in the short-period mode is of serious con- 
cern to the pilot, since the period of the oscillation can be of the 
same order of magnitude as the pilot's reaction time. The oscillation 
may therefore be difficult or even impossible for the pilot to control 
manually. Furthermore, the additional load imposed upon the airframe 
due to a rapid growth of the amplitude of a negatively damped oscillation 
makes possible the occurrence of structural failure. It is therefore of 
considerable interest to obtain an understanding of the nature of the 
short-period mode, parallel to that which has been gained of the phugoid 
mode. 


One means of viewlng the aerodynamic phenomena, occurring during the 
short-period oscillation from a fundamental standpoint is through appli- 
cation of the concept of indiclal functions. In this approach, the vår- 
tations with time of the aircraft angle of attack and angular velocity 
during the oscillation are replaced by a large number of small instanta- 
neous or step changes. The transient aerodynamic reactions to these step 
changes are termed “indicial functions," and have been calculated theo- 
retically for several classes of wings (refs. 2 to 6). By suitable 
superposition of these results (refs. 7 to 9), the aerodynamic forces and 
moments caused by the given maneuver can be studied. It will be the 
purpose of this report to make such a study for the simplified case of an 
aircraft performing single-degree-of-freedom rotary oscillations. For 
this maneuver, which corresponds to the short-period oscillation when the 
plunging velocity of the aircraft is zero, the use of simple physical 
relationships associated with the indicial function concept enables 
qualitative studies to be made of the separate effects on the aerodynamic 
forces and moments of changes in Mach number, aspect ratio, plan form, 
frequency, and thickness. Results of this investigation are then com- 
pared with the results of experiments with a group of low-aspect-ratio 
wing-body combinations. The tests were conducted in the Ames 6- by 6-foot 
supersonic wind tunnel and were similar in technique to those reported in 
reference 10. 
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NOTATION 


aspect ratio, v^ /s 
lift coefficient, lift/q58 
pitching-moment coefficient, pitching moment/qgoSc 


P-7Pg 


pressure coefficient, 


moment of inertis, slug-ft? 

free-stream Mach number, Vo/&o 

Reynolds number, based on wing mean aerodynamic chord 
wing area, including portion enclosed by body, sq ft 
flight speed, ft/sec 

speed of sound in free stream, ft/sec 

wing span, ft 

wing root chord, ft 


b/2 2 
(local chord) dy 


onm 
— 


wing mean aerodynamic chord, 


base of natural logarithms 


KT 


reduced frequency parameter, o 


local loading at plan-form surface, 
pressure lower surface - pressure upper surface 


do 
angular velocity due to pitching, radians/sec 


free-stream dynamic pressure, 3 Doo; lb/sq ft 
time, sec 
time required following an instantaneous change in angle of attack 


or angular velocity for the transient lift or moment to attain 
steady state, sec . 


— 
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x, Y, 2 rectangular coordinates 


xg. C. distance from leading edge of M. A. C. to aerodynamic center, ft 


Xo 


AXo 


Po 


Og 


distance from leading edge of M.A.C. to axis of rotation, ft 


XA. . T Xo 


angle of attack of wing center line with respect to free-stream 
direction (sketch &) 


[5-1 | 


ratio of specific heat at constant pressure to that at constant 
volume 


&irfoil-thickness ratio, maximum thickness/chord 


angle of wing center line with respect to horizontal axis 
(sketch a) 


acute angle between wing plane of symmetry and trailing edge 
(sketch p) 


free-stream density, slugs/cu ft 
angular frequency of oscillation, radians/sec 


distance traveled, measured in half M.A.C. lengths, in the time 
interval tg, AVotg/ð 


When a, d, and q are used as subscripts, a nondimensional deriv- 


ative is indicated, and this derivative is evaluated as the independent 
variable (d, d, or a) approaches zero. For example, 


ac ac 

Cm, = EJ Cm = = 

* A= o ‘po 

Vols 
& 0 
" Cm 
na i à qc 
2Vo q—0 


NACA RM AS52LOHla 9 000 2 


THEORY 


Application of Indiecial Functions to the Aerodynamic 
Theory of Unsteady Flows 


One of the most useful tools in the study of unsteady flows is the 
concept of indicial aerodynamic functions which may be defined briefly 
as the aerodynamic response of the airfoil as a function of time to an 
instantaneous change in one of the conditions determining the aerody- 
namic properties of the airfoil in a steady flow. Theoretical aerody- 
namic indicial functions were first derived by Wagner (ref. 2) for the 
two-dimensional wing in incompressible flow. More recently, these 
results have been extended by Heaslet and Lomax to cover the compressible 
case for both subsonic and supersonic speeds (ref. 4). In addition, 
theoretical indicial functions have now been obteined for both wide and 
slender triangular wings and rectangular wings, all for supersonic speeds 
(refs. 4 to 6). 


The indicial function derives its usefulness primarily through the 
ease with which it lends itself to the powerful and well-established 
methods of the operational calculus (refs.7 to 9). With the use of these 
methods, the aerodynamic forces and moments caused by arbitrary motions 
of the airframe can be studied from a fundamental standpoint. Because of 
the wide range of applicability of this means of approach in unsteady flow 
analyses, a considerable portion of the succeeding discussion is devoted 
to the fundamentals involved. 


Definition of coordinate system.- In the succeeding analysis the 
stebility system of axes is used. The origin of the coordinate system 


is placed in the airfoil so that the y axis which is perpendicular to 
the vertical plane of symmetry is coincident with the axis of rotation of 
the airfoil; the positive branch of the x axis is pointed in the direc- 
tion of flight; and the z axis lies in the vertical plane of symmetry, 
positive downward. The angle of attack a is measured as the angle 
between the chord plane of the airfoil and the xy plane, and is shown 
as positive in sketch (a). The angle of pitch @ is the angle between 
the chord plane of the airfoil 
and the horizontal plane (an 
arbitrary reference) and is also 
shown positive in sketch (a). 
Forces are measured as positive 
upward, whereas pitching moments 
are positive when tending to 
increase the angle of pitch in 
the positive direction. When 
the airspeed vo is constant 
which corresponds to the 


flight path \ 


V 


Sketch (a) 


— . 
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condition under study, the translatory and angular motions of the air- 
foil with respect to any system of coordinates are defined if the time 
histories of the angle of attack a and the angle of pitch 0 and 
their derivatives are known. For purposes of clarity, two different 
harmonic motions of the aircraft are shown, illustrating the difference 
between a flight path which involves a constant angle of attack and a 
varying angle of pitch (sketch (b)) and one which involves a constant 
angle of pitch and a varying angle of attack (sketch (c)). 


2 , 
flight path flight path 
angle of pitch = 8 angle of pitch =0 
angle of attack = Q angle of attack & 
Sketch (b) Sketch (c) 


Now consider the case of a wing executing harmonic rotary oscills- 
tions about the y axis while the origin of the coordinate system 
traverses a level path at constant velocity vo. This case corresponds 
to that of a wind-tunnel model mounted to permit single-degree-of- 
freedom rotary oscillations, or to the short-period mode of an aircraft 
in flight when the plunging velocity of the center of gravity is zero. 
Here a and O are equal, so that the maneuver is defined by one vari- 
able, the time history of either a or O. Let the angle of attack 
be a and the angular velocity be q (q = de/dt = da/dt). At any 
instant, the normal velocity at any point on the airfoil surface is com- 
posed of. two parts, one due to the instantaneous angle of attack  aVo, 
the other due to the angular velocity at the same instant -x (see 
sketch (d)). These are two of the instantaneous boundary conditions of 


the unsteady flow. 


normal velocity Va 
due fo angle of attack a 


normal velocity -gx due fo 
angular velocity g 


Sketch (d) 
Seer 


Ll 


* 
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Solutions for the aerodynamic forces and moments which correspond 
to these boundary conditions may be derived by & number of methods 
involving various degrees of approximation. In succeeding sections, 
the use of the concept of indicial functions and the principle of super- 
position for this purpose will be illustrated &nd compared with other 
current widely used methods. 


Concept of indicial functions.- In order to illustrate this concept, 
assume that the airfoil under consideration has been flying a level path 
at zero angle of attack. At some time, which is designated time zero, 
the wing is caused to attain simultaneously a constant angle of attack a 
and angular velocity q. The normal velocity of the flow next to the 
surface of the airfoil therefore changes discontinuously from zero to & 
pattern that is constant with time and identical in shape to the pattern 
shown previously in sketch (d). The lift and pitching moment that 
result are of a transient character and attain their steady-state values 
corresponding to these new boundary conditions only after a significant 
interval of time has passed. It should be noted there exists an essen- 
tial difference between the length of this time interval at subsonic and 
supersonic speeds. At supersonic speeds, the vorticity shed into the 
airfoil wake cannot influence the flow about the airfoil but at subsonic 
speeds this influence exists for all time. The result is that the lift 
and moment reach steady-state values in a finite time at supersonic 
speeds but approach these values asymptotically at subsonic speeds. In 
either case, however, the time responses in lift and moment to the step 
changes in a and q are termed indicial functions. Sketch (e) illus- 
trates typical subsonic and supersonic indicial lift responses to a step 
change in the angle of attack. 


C, M>! 


C, Clg) 


Ip } 


Sketch (e) 
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It is obvious that the time history of the wing motion during a | 
short-period oscillation may be broken down into an infinite number of 
infinitesimally small step changes in the angle of attack and step 
changes in the angular velocity. The summation of the indicial lift and 
moment for these steps then yields the total lift and moment at any pre- 
seribed time. In sketch (f), the mechanics of the procedure are illus- 
trated for an arbitrary angle-of-attack variation. Here, the given 
angie-of-attack variation is 
replaced by a number of small 
step changes. Within each step 
the corresponding response in 
Lift is shown plotted for con- 
venience. It is then apparent 
that the total lift at time t 
is equal to the sum of the 
increments of lift in each step 
at time t. As indicated by 
the leaders, however, it is 
clear that the increments of 
lift for the various steps at 
time t are equivalent to 
increments in the first step 
at time t-t,. Alterna- 
tively, then, the total lift 
at time t can be written ag! 


Sketch (f) 
t 

Cy(t) = Cx4 (*)o(0) + J Cra C t- ti) 185 (1) Ab (1) 
O 


After & transformation of variables, t - t, - T, and letting the incre- 
ment of time approach zero, equation (1) caen be rewritten in a form of 
Duhamel's integral (see, e.g., ref. 9) 


Oplt) = f° ole aT (2) 


+Here, and in the remainder of this report, the use of parentheses is 
reserved solely for the indication of functional dependence. Thus, 
for example, in equation (1) the term Cra(t-ti) is interpreted as 


the value of Che at time t - tj, whereas AE (ta) is the value 
of Mut at time ti. All other enclosures indicate algebraic 
expressions in the usual sense. - Í 
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A similar procedure is carried out for the angular velocity 
variation, whereupon the total Lift coefficient at the prescribed 
time t becomes 


x 


C(t) = F Cr (T)a(t- ar + — a Cn (r) ay „ a(*-T)aT (3) 


It should be pointed out that in this form equation (3) is applicable 
to the analysis of arbitrary motions, the only restriction being that 
the flight speed is constant. In the following sections, however, the 
application of equation (3) is restricted to harmonic motions having a 
single degree of freedom. The reasons for this restriction are two- 
fold: first, the motions of a statically stable aircraft in response to 
a disturbance are most generally of a harmonic nature; and second, such 
a restriction permits an assessment of the influence of the time rate 
of the airfoil motions on the aerodynamic forces and moments. 


Application of indicial functions to harmonic pitching oscillations - 
Consider first a pure sinusoidal pitching oscillation, the angle of 


attack being zero throughout the motion. The flight path for such a 
motion has been illustrated in sketch (b). In this case, the angle of 
pitch is given by 


alt = gie a" 


where O09 is the maximum amplitude of oscillation and w is the 
angular frequency. The angular velocity is, of course, q(t) = 6 = 
108,6 eit = iw9(t). Inserting the value for q(t) in equation (3) and 
per tormine tne indicated operations, there results 


Cr t) = = uz o(t) L Cr Cr)e7 rar To—— = 


Vo 7, 901409) (4) 
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Note in sketch (g) that Cra C7) is equal to Cra (t) FZ (r), and that for 
subsonic speeds Fs(tT) approaches zero as 7 approaches t, 


R>) 


G (r) Cog lt) C (02) 
= 


Sketch (g) 
Replacing Cra(7) in equation (4) by this equality, 


C(t) 100 wrt t -iuT 
"ME m Cr C8) E [ Fo(Tr)e ^" àv (5) 


For subsonic speeds, let t approach infinity. With this sub- 
gtitution, equation (5) thereby represents the lift coefficient due to 
the harmonic pitching motion after the transient loading subsequent to 
the start of the motion has reached steady state. Then separating 
equation (5) into components in-phase (real part) and out-of-phase 
(imaginary part) with @, there is obtained 


* 2 A n Fe (r) cos wrdt 
(6) 
Cy, B on 
75897 = Cr, (e) - W 5 PS(T) sin wrdt 


Introduce the nondimensional parameters, 


2 

S — T number of half M.A.C. lengths traveled in time T 
we - 

k = — reduced frequency 
2Vo 
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In terms of these parameters, equation (6) becomes, for Mo 1, 


Sb K f Fe (ꝙ) cos kd 

9 o 
. (T) 
L eo 

— — (Cl — k F in kod 

qé/2V, 14l) js. bah) Bin Kpg 


At supersonic speeds, equations (T) may be simplified somewhat 
since the build-up in lift is completed in a finite number of half 
M.A.C. lengths of travel og. In equations (7), therefore, the upper 
limits of the integrals may be replaced by og, since beyond that point 
FS(Q) is zero. 


c 
M k= f * Fa (9) cos Kd ꝙ 
8 0 
P 2% (8) 
Mo > 1 
hor (oa) -k[ "FrQ(9) sin pap 
q8/2Vs 9 T O 


Thus, it appears from equations (7) and (8) that there are both in-phase 
and out-of-phase lift forces associated with the harmonic pitching oscil- 
lation. Notice, however, in equations (8) that if the trigonometric 
terms are expanded and the reduced frequency is required to be small 
compared to unity (corresponding to the frequencies encountered in 
dynamic stability work) terms containing second and higher powers of k 
will be very small compared to first-order terms. Thus, for slow fre- 
quencies, the only force of consequence during the pitching oscillation 
is the first order in frequency out-of-phase lift force, dqẽ / 2 Cro (oa) ^ 


2 This quantity is, of course, the same lift force in phase with the 
pitching velocity which would occur alone had the wing been executing 
& steady turn (q constant). One of the chief advantages of the indi- 
cial response method, at least for supersonic speeds, is the ease with 
which the relative importance of the various terms contributing to the 
total lift and moment can be assessed and the sources of the important 
contributions identified. 


— — 
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The phase relationships for the har- 


q monic pitching oscillation are indi- 
ge C, (o) cated in sketch (h). It is evident 
eu ^? N from sketch (h) that the total lift 


leads the angle of pitch by nearly 90°. 


C (total) Application of indicial functions 
to harmonic plunging oscillations.- 


Next, consider a purely sinusoidal 
variation of the angle of attack, the 
55 — MÀ angle of pitch being zero throughout 
y. TT 7 ihe motion. The flight path for this 
- ge "d sinkgdg ok Fluß cost og motion has been illustrated in 
£V, "o á sketch (c). Here, a equals apet, 
where, as previously, ao and w are 
the maximum amplituđe and angular fre- 
quency, respectively. Applying equa- 


Sketch (h) tion (3) again, 
t ee 
Cult) = ao dc J, Cial?) mena 
x i ea Or, (7) e" Tar + acr, (t) (9) 
o 


Now, as in the previous example, let Cr, (T) = Cr C9) - F (T) so that 
equation (9) becomes 


Cr(t) = a(t)cr (t) - ih, t) p Fi(7)e719" ar (10) 


Again, introduce the nondimensional parameters Ọ and k, separate 
equation (10) into its real and imaginary parts, and let ꝙ approach 
infinity for subsonic speeds and og for supersonic speeds. There 
results | 


CT, 
G 


X 
= On - x f Fi (S) sin kpdg 


o^ 
| 
3 
S 
A 
I 


(11) 


> 
it 
Q 

© 

*e 

ce 
V 
H 


Cr * | 
87e = a ij F, (9) cos Kk 
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Notice in equations (11) for supersonic speeds that when the trigono- 
metric terms are expanded for the slow frequency case, as was done in 
the previous example, there appears an in-phase term of zero order 

in k, aCp,(0g), and an out-of-phase term of first order in k, 


. C 
- aL f "E()ap. These, then, are the principal contributions to the 
O 


lift forces for the slow-frequency angle-of-attack variation. The phase 
relationships for this motion are shown graphically in sketch (i). For 
this case, it is evident that the total-lift force can lag behind the 
angle of attack because of 
the negative out-of-phase 
contribution, & 
M % 
2:6 p Fa (S) cos kpdo - ak cen toes «C, (a) 


2b "o A rd s 


Application of indicial 


functions to harmonic rotary 

oscillations.- Finally, con- 

sider the case of harmonic oa 

rotary oscillations. Here, — &c Ti ipleos kde 
as previously mentioned, the S 

normal velocity over the 

wing surface is composed of | 

contributions from both the 

angular velocity and the Sketch (i) 

instantaneous angie of 

attack, so that the complete 

expression in equation (3) must be employed to obtain the totel Iff. 
However, for single-degree-of-freedom rotary oscillations, œ equals O0 
and d equals q, so that in this case the separate expressions given 
for the harmonically pitching wing (eqs. (7) and (8)) and the harmoni- 
cally plunging wing (eq. (11)) can be combined to give the total lift 
for a wing executing harmonic rotary oscillations. Then, adding the 
results of equations (7), (8), and (11), the in-phase and out-of-phase 
components of the total lift become 


G, (total) 


= Cn) -k 18 i Fi (SH) sin kpdg + k? f 5 Fo(®) cos kQpdg (12a) 
Xx, MSL 
NO, M92 1 
Cr, P À 
"774 = Cra) -k 5 Fe (S) sin k9aQ - A Fi(9) cos k9ao (12b) 
O 


— ̃— 
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The phase relationships for the rotary oscillation may, of course, also 
be obtained by directly adding the results given in sketches (h) and (i). 
The result of this addition is shown in sketch (j). 


q A 
2 
$ a Kg (Blcos bes 


A 
8 nkd d. 
z f E Via 7i 


A 
- g [Fi igicoskod¢ 
2W? 


A 
ela (a) E Helen tede) 


o 


Sketch (3) 


It will be noted in sketch (3) that the total-lift force can either lag 
behind or lead the angle of attack, depending on the relative magnitudes 
of the three terms comprising the out-of-phase lift. The total lift is 
shown lagging behind the angle of attack in sketch (j), which situation, 
for axis positions ahead of the point of concentration of the total lift, 
gives rise to the possibility of the development of negatively damped 


rotary oscillations. 


Again, the complete frequency-dependent equations for the total 
lift of & wing in supersonic flow due to the rotary oscillation 
(eqs. (12), X = da) may be reduced to first order in k for the slow- 
frequency case in the same manner as was described in the two previous 


examples to give 


CL 

— = Cz,(0a) 
: (13) 
L n | r9& 

"HEP Clg (oa) J Fi(q)ào 


3X 
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For all three examples, the same procedure may, of course, be used 
to obtain the pitching-moment coefficient. Only the pitching-moment 
equations for the rotary oscillation are presented below, since the 
correspondence between the lift and moment equations 1s obvious. For 
the rotary oscillation case, then, 


. 


A A 
Cm . Cm 0M) = * 6 FS(9) sin kPdP + E Fa (S) cos kPa? 


= <i 
ROA, Mg» 1 


= = Cm, (X) -K [ F. (S) sin kpdp - D Fe (P) cos kpap 


where, as previously, 


Fs(9) = Cm) - Cn (9) 


and 


F4(®) = Cng(d) - eng () 
Again, reducing equations (Ih) for supersonic speeds (A=0,) to first order 
in frequency, there results 


zm o. Cm KC O 
(15) 


C o 
m & p 


qE/ANS Cmg( oa) - L s (9) ap 


The complete frequency-dependent equations for the lift and pitching- 
moment coefficients for the rotary-oscillation case (eqs. (12) and (14)) 
describe completely the aerodynamic forces and moments resulting from 
the single-degree-of-freedom pitching mode. For the purposes of the 
present discussion, however, it is sufficient to limit consideration to 
the simpler first order in frequency results of equations (13) and (15). 
The significance of the effect of the higher-order terms on the out-of- 
phase pitching-moment will be examined in a later section of this report. 


With regard to the first-order results for the supersonic lift and 
pitching-moment coefficients, it is instructive to note that the 
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g 0] 
quantities f > r (S) dy and L > Fg(9)àQ in equations (13) and (15) 


are d gecnetsicuiis by the &reas of the shaded portions of 
sketch (k). 


Sketch (k) 


The manner in which these areas are affected by variations in Mach 
number, aspect ratio, plan-form shape, and thickness will be used as a 
guide in later sections of this report to determine the significance of 
these parameters. 


In the foregoing discussion, no mention has been made of reducing 
the complete equations for the lift and moment coefficients at subsonic 
speed (eqs. (12) and (14), X=) to first order in frequency as wes done 
for the equations noted as applying at supersonic speed. It is evident 
that if the same procedure had been applied for subsonic speeds, the 
area, corresponding to the term f Fi (S) dꝙ can either be finite or can 
become infinitely large, depending on the manner in which the indicial _ 
lift function OCr,(9) approaches its steady-state asymptote as 9 >o. 
In the latter case, there exists the interesting anamoly of an infinite 
out-of-phase lift force as the frequency approaches zero. As can be 
seen from the results of reference 4, such will be the case for the two- 
dimensional wing. This result as the frequency approaches zero is not 
peculiar to the indicial analysis alone, but has been pointed out by a 
number of authors using different approaches. As indicated by Miles 
in reference 11, however, the anamolous result can be considered to be 
& consequence of assuming a two-dimensional flow, and there is reason 
to believe that the difficulty as. the frequency approaches zero will not 
exist for finite-span wings. 


Ás a been W F the use of the characteristic 


areas f Fi (S) dꝙ and Fo Fg (O) do will be shown to be of considerable 
value in estimating the damping-in-pitch characteristics of wings at 


supersonic speeds. For the two-dimensional wing at subsonic speeds, the 
singularity as k~O prevents the use of such a simplified approach 
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without further study. However, rather than return to the use of the 
full frequency-dependent equations ((12) and (14)), the reduction of 

the eguations for the out-of-phase lift and moment to first order in 
frequency will be made in such a manner as to preserve the significance 
of these areas. To accomplish this end, equation (12b) is reconsidered. 
It is evident that the first integral in equation (12b) may be discarded, 
since its contribution to the out-of-phase lift is at least of second 
order in frequency. The second integral is divided into two parts: 


eo co 

f Fi(9)ceos Kd = ig Fi(9)eos R + f Fa (o) co kPdP (16) 

2 | Q4 
where 9, is chosen such that FI (O) is close to zero. The first 
integral in equation (16), being bounded, then causes no difficulty. 
Expanding the trigonometric aor and retaining only the first term in 
the expansion, there results ip T1 P. (S) d, which is tze characteristic 
area out to the point Qi. Now “For large values of 9, Fi (O) is approx- 
imated in reference 4 by 


Eat?) = TE - p` um 


where the values of h, a, and b are dependent on Mach number, and are 
given for Mo = O, 0.5, and 0.8 in reference 4. Inserting this quantity 
in the second integral in equation (16), we have 


eo 9€ cos kP cos K 
es. Kd = u Lato do + ub f. Teo a9 (17) 


Performing the indicated integrations in equation (17), there is obtained 
a term, ban, from the second integral, and a term, -uüCi4k[sacQ, ] 


from the first integral. For small values of the argument, the cosine 
integral is &pproximated by (see ref. 12), 


Ci L [8-04 1} e e [8-94 1} 


k[849,] — 0 
where y is Euler's constant, 1.78107. 


Then, through the first order in frequency, 


Ts 


= Chg (») - Fi (S) dp + G(9, „K) (18a) 


2 
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where 


G(9,,k) = [um piteeoi1 - Bo | 


a+ Q4 


The out-of-phase moment for pitching &bout the leading edge follows from 
the above development, with the added result of reference 4 that for 
large values of Pp, 


F; 
Fs (ꝙ) z - ne 
Then, through the first order, 
C 9 
ar, nd) - fj reg - F Gd, (185 


Thus, after fixing 031, choosing a (small) value of k, and computing 

G(?, „k), the finite areas corresponding to the terms f Pı F, (P)ae 
0 

and - f 1 FQ(9)d? can be assessed in the same manner as will be done 


for the supersonic case. The advantages of such a procedure will be 
evident later. 


Correspondence between indicial lift and moment analysis and other 
methods. Before proceeding further with applications of the indicial 
response method, it is appropriate to discuss the relationship of this 
approach to other widely used methods. 


Following the fundamental papers of Bryan and Routh, which intro- 
duced the basic differential equations of motion of rigid bodies and 
their stability criteria, the historical development of the theory of 
longitudinal motions of an aircraft evolved separately in two fields of 
research: dynamic stability and flutter. Workers in the dynamic stabil- 
ity field soon found that the longitudinal oscillations of a rigid air- 
craft in flight were generally of small reduced frequency. On this 
basis, the constants due to the aerodynamic properties of the airframe 
which appear in the differential equations of motion were considered to 
be independent of frequency. As a first approach to the problem of 
obtaining the necessary aerodynamic coefficients analytically, the 
instantaneous normal velocity distribution at the surface of the airfoil 
was assumed to be constant with time. The aerodynamic forces and moments 
arising from the fixed boundary conditions were then calculated using 
steady-flow theory. Later, this assumption was realized to be an over- 
simplification for the case of wing-tail combinations and an additional 
term correct to the first order in frequency was added which accounted 
for the lag in the tail pitching moment caused by the time required for 
the vorticity discharge from the wing to reach the tail (see ref. 13). 


— 
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Since at low speeds the pitching moment of the teil far outweighs all 
other contributions, the results from steady-flow theory together with 
the term accounting for the vorticity lag satisfactorily predicted the 
dynamic longitudinal motions of wing-tail combinations, and it was con- 
eluded that the major aerodynamic effects had been accounted for. In 
recent years, however, numerous authors (in particular, Miles; see, e.g., 
ref. 11) have pointed out that the above-mentioned theory overlooks 
important contributions to the aerodynamic forces and moments which, 
though still within the first order in frequency approximation, arise 
from time-dependent boundary conditions and must be calculated from 
unsteady-flow theory. It has been shown by these authors that with 
proper inclusion in the equations of motion of these coefficients, the 
deterioration of damping in the short-period mode actually occurring for 
aircraft flying at speeds near the speed of sound can be successfully 
predicted. The consequences of the assumptions involved in the classical 
dynamic stability theory will be more evident from a brief review of the 
equation of motion and boundary conditions for the single-degree-of- 
freedom rotary oscillations of a rigid wing flying at constant supersonic 
speed. At the very outset, the assumption is generally made that the 
aerodynamic reactions to the motion of the airframe depend only on the 
angular position and angular velocity and not upon angular accelerations 
or higher time derivatives. The equation of motion for the change in 
pitching moment following a displacement from an equilibrium position is 
then written in the form of a power series: 


Is EIS zm veiba E, [Ea] [As], 
9085 3 [hc 72701 x 3 [ac rm dps ov, Ld 
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It should be remembered that for the rotary-oscillation case, the air- 
foil is subjected to changes in both angle of attack a and angular 
velocity q, and that these motions produce normal velocity patterns at 
the airfoil surface which are different in character. Thus, although 
for the single-degree-of-freedom case, d and q are equal, nevertheless 
their separate effects must be considered and it is therefore 

8 [3e/2v] 3 [G6 /2V4 ] 
Next, if it is assumed that the moments are linearly dependent on their 
respective variables, the higher-order terms in equation (19) may be dis- 
carded and the remaining partial derivatives considered as constants for 
the given wing. There remains, therefore, & linear second-order system 
with constant coefficients. In order to calculate the coefficients 


— 


necessary to include both in equation (19). 
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(termed stability derivatives) theoretically it became necessary, for 

lack of more refined theoretical methods, to assume that the instanta- 

neous normal velocity of the flow at the surface of the wing was fixed 

with respect to time. Thus, the partial derivative PE ris NES could 
SLꝗ / vo. 

be calculated as the pitching moment due to a constant pitching rate, 

that is, COng( 98) » while the derivative OC,/da becomes the pitching 


moment due to a constant angle of attack, that is, Cm (08). As a con- 


sequence of fixing the normal velocity pattern in time, however, it was 


necessary to assume that the derivative 2 was zero. There 


39 [&c/2vo ] 
was therefore no possibility for this theory to predict the occurrence 
of dynamic instability for a wing alone, since the only damping term 


remaining is Cmꝗ (o), which is always stabilizing. When the restriction 
of constant normal velocity with time ig lifted, however, the assumption 
is then made that the stability derivatives in equation (19) may be caleu- 
dated separately by fixing each of the independent variables a, &, 

and q in turn with respect to time. The derivatives Cmq oa) 
—— (or Cp. 
S LGS / ev 0 me) 


can now be included and calculated as the pitching moment due to a con- 
stant vertical acceleration, aV,. It should be emphasized that while 
Cm (oa) and Cma(c may be calculated from steady-flow theory by virtue 


of the assumed invariance with time of the normal velocity pattern, Cmó, 
must be calculated from unsteady-flow theory since for constant &. 

the angle of attack varies linearly with time, as does the normal veloc- 
ity of the flow &t the surface. | | 


and Cy (oa) thus remain unchanged, but the derivative 


It is clear that since the stability derivatives in equation (19) are 
assumed to be independent of the frequency, the result for the aerody- 
namic pitching-moment coefficient is thereby limited to one that is 
correct only to the first order in frequency. For the single-degree-of- 
freedom case, then, the in-phase and out-of-phase components of the 
total aerodynamic pitching-moment coefficient correct to the first order 
in frequency, become 


=n = Omer Dg) 
(20) 
Cm 


MJER = Cng(ca) Cm (O 


By comparison with the first order in frequency result from the indicial 
response analysis (eq. (15)), it is evident that the two results are 
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G 
identical if the quantity - f = Fg (O) dp can be shown to be equivalent 
O 


to Cys(og). To show this equivalence, consider a wing, initially in 
level steady flight, which is suddenly forced down with constant vertical 
acceleration WO. As seen in sketch (1), the angle-of-attack variation 
in this case is a = åt, where d. is a constant. Then applying the 
counterpart of equation (3) for the 

pitching moment 


Cy(t) = = 05 Cm, (7) d t- T) dt 


Inserting a(t-T) = &lt-T], and per- 
forming the indicated differentiation, 


C(t) = &. p cm T) dr 


Now replace Cm(T) by Cm4(t)-Fa(7T) o (0) Wt) / 
and let t be greater than tg. 
Then Sketch (1) 


Cm(t) = a(t)Cm (ta) - & in Fg (T) d 


and nondimensionalizing, by replacing t and tg by co / Vo and cg / Vo- 
we have 


Cm(9) = a(9)Cn, (ca) - — ~ Fs (S) a5 P 2 0 (21) 
O 


Thus, the pitching moment proportional to the constant vertical acceler- 
ation parameter G€/2V,, which is synonomous with the definition of the 
stability derivative Cys, is found to be equivalent to the pitching- 
moment contribution due to d for the first order in frequency-rotary- 
oscillation case. Therefore, the results of the indicial response 
method, when reduced to the first order in frequency for supersonic 
speeds, are identical to the results from the familiar first-order theory 
used in dynamic stability work.* 


Spy the same procedure, the stability derivative CIS (on) can be shown 


to be equivalent to the term - [5 F, (9) dp. 


“Notice in equation (21) that if o, is replaced by infinity, the results 
apply to subsonic speeds. For the two-dimensional wing, the analogy 


e 
between Cm; end - 3 Fg (ꝰ) dd then gives only the previously mentioned 
singular i ty at infinity as Kk — O. If the area corresponding to 


f F,(P)d® were finite, however, the analogy would be equally useful 


for subsonic as well as supersonic speeds. 
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Workers in the field of flutter, who were concerned with frequencies 
many times those encountered in dynamic stability analyses, required 
theoretícal information showing the behavior of the forces and moments 
as affected by the frequency of oscillation, and therefore discarded the 
first-order theory for more precise methods of analysis. One of the 
most useful of these has been the "oscillating potential” theory, which 
is based on periodic solutions of the time-dependent linearized equation 
of compressible potential flow. The in-phase and out-of-phase lift and 
moments are thereby determined, generally as functions of powers of the 
reduced frequency, aspect ratio, Mach number, and position of the axis 
of rotation. The application of this method, which developed primarily 
ag a result of Theodorsen's work for incompressible flow (ref. 14) has 
recently produced a number of useful papers covering a wide variety of 
wings at supersonic speeds (see, e.g., refs. 15 to 19}. It has been 
shown by a number of authors, in particular, Garrick, in reference 7, 
that through the use of superposition methods the results for the aero- 
dynamic coefficients obtained from the oscillating potential theory are 
wholly compatible with those of the indicial response method for har- 
monic motions (eqs. (12) and (14)). 


Thus the indicial response method embraces both the first-order 
theory of dynamic stability and the oscillating potential theory, and, 
in effect, bridges the gap between the fields of dynamic stability and 
flutter. Furthermore, the indicial response method overcomes the dis- 
advantages of the two methods described above since, unlike the first- 
order theory, the effects of frequency on the aerodynamic coefficients 
can be determined, and, unlike the oscillating potential theory, the 
method can be applied easily to the study of arbitrary motions. Finally, 
the indicial response method represents a fundamental approach to the 
problem of unsteady flow, and affords valuable insight into the physical 
nature of the aerodynamic phenomena taking place. 


Physical concepts relating to the indicial loading.- It has been 
shown that for even small frequencies, the pitching moment of an airfoil 
in harmonic rotary motion can lag behind the angle of attack of the air- 
foil. The magnitude of the lag depends on the character of the indicial 
response to a step change in angle of attack. It is therefore of 
interest at this time to re-examine the physical nature of the flow 


that contributes this lag. 


Consider first the lift and moment at the instant the angle of 
attack changes, assuming that previous to time zero the wing has been 
flying a level path at zero angle of attack. At t = O, the wing begins 
to sink, without pitching, with constant downward velocity q while 
maintaining its forward velocity. The angle of attack therefore changes 
discontinuously from zero to a constant a. At the same instant, the 
step change in angle of attack causes the emission of a compression wave 
from each point on the lower surface of the wing and expansion waves 


+X 
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from points on the upper surface. In the infinitesimal time during which 
the starting action occurs, each section of the wing experiences the 

same impulsive force, and by equating the impulse to the momentum trans- 
mitted to the mass of fluid affected by the starting waves, the starting 
lift coefficient can easily be derived as 10% ' (see ref. 3). During 
the infinitesimal starting time, the pressure disturbances from the edges 
of the airfoil, propagated at the speed of sound, travel an insignificant 
distance and do not influence the remainder of the airfoil. The lift 
coefficient is therefore independent of the wing plan form. This remark- 
&bly simple result for the starting lift coefficient, which is valid for 
both subsonic and supersonic speeds, is thus dependent solely on the 
flight Mach number. The starting pitching moment follows directly from 
the above result, since by virtue of the uniformity of loading the aero- 
dynamic center is located at the wing centroid of area. 


For values of time greater than zero, however, the situation differs 
radically for the supersonic and subsonic speed ranges. Consider first 
the supersonic case. As time passes, the spherical sound waves emitted 
at t =0 grow in size with radius apt. The wing, however, is moving 
forward at a faster rate than the rate of growth of the starting sound 
waves and thus begins to emerge from the influence of these waves. This 


is shown schematically in sketch (m). leading edge at t=t 
= fg 


Starting sound waves 
from leading edge 


a? f» ]j, 
envelope of s 


starting waves 
from leading 
edge of tz 


leading edge a! t= 
leading edge a? t=0 


FFF 


Sketch (m) 


At t = O, the starting waves just cover the wing and the loading is 
uniform as described previously. At t = t,, the starting waves have 
grown in radius and the wing has begun to emerge from their influence. 
On that portion of the wing which has emerged, region (1) in sketch (m), 
the loading has already reached its steady-state value. Notice that in 
this region the characteristic tip Mach cone has already formed. On the 
portion of the wing uninfluenced by the starting waves from the edges, 
region (3) in sketch (m), the loading is still uniform as at t - O. 

In regions (2) and (4), the loading is influenced by the starting waves 


. 
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from the leading and side edges, and in these regions is thus different 
from the loading in either regions (1) or (3). As time increases still 
further, the uniform starting load quickly disappears as the sound waves 
from the leading edge grow in size and as the wing moves forward. | 
Finally, at time ta, the envelope of the starting waves from the lead- 
ing edge is coincident with the trailing edge of the wing, and the 
steady-state loading corresponding to the new angle of attack a has 
been completely established over the wing. 7 
The above relationships can be shown more clearly for the entire 
time interval zero to tą for a two-dimensional wing by plotting as a 
function of time the position of tbe wing leading and trailing edges 
and position of the envelopes of the sound waves which emanate from the 
leading and trailing edges at t - O. Such a plot is shown in sketch (n). 


frace of 
leading edge 


W 7 trace of starting E 


frace of sound waves from 
frailing edge | /eading edge 
7 
Sketch (n) 


It is clear that at t =t, the regions of the wing (1), (2), 
and (3) correspond to the same regions at t = tı for the wing shown 
in sketch (m). For t =0 and in region (3) the loading is uniform and 
is given by lha/Mg. For ttg and in region (1) the wing has out- 
stripped the starting waves from the leading edge and has attained its 
steady-state loading. For t<tg the chordwise loading is composed of 
combinations of the loading in each of the three regions as shown in 
sketch (0).9 


SThe reader will note the similarity between sketch (n) and sketches 
depicting the boundary conditions for three-dimensional wings in 
steady supersonic flow. Many researchers have pointed out the analogy 
and it has been used to calculate the pressure over a wing impulsively 
starting from rest (refs. 3 and 4). 
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Since the loading on the 

wing attains its final 
steady-state distribution 

at precisely the time when 
the wing has emerged entirely 
from the influence of the 
starting sound waves from 

the leading edge (or apex), 
the time to reach steady 
state may be easily calculated 
for any type of wing by means 
of the geometric relationships 
shown in sketch (p). 


Sketch (o) 


leading edge 
at f= fa 


trace af t= fa 
of starting sound wave 
from opex 


leading edge 
af 0 


Sketch (p) 


e 


26 CONT Pa. | NACA RM A52L04a 


lt may be easily verified that ta, the time required for the wing to 
attain its steady-state loading, is given by the following relationship: 


See Inc WM y< (22) 
VolMo-ese v] 2 


Notice that for wings having straight or swept-forward trailing edges 
and straight or swept-back leading edges, equation (22) reduces to 


= — Mo _ 2 
ta vol 11 (23) 


In terms of the number of half M. A. C. lengths of travel, it thus appears 


that for wings having straight or swept-forward trailing edges and 
straight or swept-back leading edges, 


(24) 


whereas for wings having swept-back trailing edges and e or 
swept-back leading edges, 


== —— 1 EE NO AN UR RE > — 
Og = -= M- y Mo 2csc V + v cos V 


— — : (25) 
Ca = Mo = Mo + 1 -e [M] ; Mo Sese v + — cos v 
c 


Mo? - 1 


The second of equations (25) applies to that range of Mach numbers for 
which the trace of the starting sound wave from the apex is not tangent 
to the trailing edge at t = tg. 


Now consider the subsonic case. Here the situation is more compli- 
cated in that since the starting sound waves travel faster than the 
wing, the wing never escapes their influence. Furthermore, the vortic- 
ity shed by the wing at t O can also influence the loading on the 
wing since the disturbances created by the shed vorticity travel forward 
at a faster rate than the wing. For these reasons, the indicial loading 
at subsonic speeds approaches its steady-state distribution asymptoti- 
cally with time. The situation for subsonic speeds will be more clearly 
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evident from examination of sketch (q), which shows the relationship of 
the traces versus time of the leading and trailing edges of a two- 
dimensional wing flying at a subsonic Mach number to the traces of the 
envelopes of the starting sound waves. 


no 
N 
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Sketch (q) 
Notice in sketch (q) that the starting 
sound waves intersect the edges of the wing 
and that each intersection causes a new 
sound wave to be emitted, which in turn will 
intersect an edge. Furthermore, notice that 
the vorticity shed from the trailing edge 
at time zero can influence that portion of 
the wing behind the sound wave trace 
. Labeled. a-b. The influence of each succes- 
sive sound wave reflection, however, is 
weaker than the last, and as the wing moves 
away from the starting vortices their 
influence diminishes, so that at t =o the 
loading on the wing attains its steady-state 
distribution. The variation of the chord- 
wise loading with time for the two- 
dimensional wing flying at a subsonic Mach 
number is shown in sketch (r). Notice that 
for eo, the chordwise loading is markedly 
different from the loading at supersonic 
speeds (sketch (o)). However, for t =0 


Sketch (r) 


aE pose 
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and in the region corresponding to region (3) of sketch (q), the loading 
is uniform and equal to ͥ , as in the supersonic case. 


Damping in Pitch of Low-Aspect-Ratio Wings 


Previously (ed. (20)), it was shown that for single-degree-of- 
freedom, slow-frequency, pitching oscillations of a wing, the principal 
parameter contributing to the damping of the motion is the damping coef- 
ficient Cy, + Cn. This result, however, is not directly applicable to 
the case of an aircraft in flight, since generally additional damping is 
provided by virtue of the fact that the aircraft experiences harmonic 
vertical translatory oscillations as well as the rotary oscillations. 

For the slow-frequency case, however, the effect of the second degree of 
freedom on the damping can be determined by use of the relationships 
given in reference 10. It is further shown in reference 10 that although 
the effect of the translatory oscillation is usually to increase the 
total damping, nevertheless, the parameter of primary importance in 
determining the magnitude and duration of the motion remains the damping- 
in-piteh coefficient Cinq + Cms. The remainder of this section is 
therefore concerned with' a stad of the effect on this parameter of 
certain important variables. In particular, the effect of the position 
of the center of gravity, and the effects of Mach number, aspect ratio, 
plan-form shape, frequency, and thickness will be examined, principally 
by inspection of the indicial lift and moment responses to a change in 
angle of attack. 


Effect of static margin.- From the previous discussion, it will be 
remembered that at supersonic speeds the stability derivatives CIA (O 
and Cm oa) were shown to be equivalent to the indicial lift and moment 


0 0 
expressions, - f > F (p)ap and - f e Fa (o) a0. Furthermore, it was shown 


that these were the contributions which could cause the total lift and 
moment during the short-period oscillation to lag behind the angle of 
attack. Hence, by inspection of equation (15), it is evident that since 
Cg (Oa) is always stabilizing, when Cy,(cqa)is negative (corresponding 
to & statically stable condition) the possibility of dynamic instability 
in the form of negatively damped rotary oscillations arises when 

Cuz (og), the shaded area in sketch (k-ii), is larger than Cmnq () 

Now since the normal velocity at the surface of the wing due to the 
instantaneous angle of attack is constant over the wing, the lift 
derivative Cri; which arises from this boundary condition is independ- 
ent of the axis position and Cys will therefore vary linearly with 
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axis position. This variation is illustrated in sketch (s), where Axo 
represents the distance of the center of gravity from the aerodynamic 
center, measured positive forward of the aerodynamic center. 
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Sketch (s) 


The parameter Cy,,, on the other hand, is a direct function of the 

axis position, since it arises from a normal velocity distribution that 
varies directly as the distance from the axis. The moment coeffi- 

clent C will therefore vary as the second power of the axis posi- 
tion, and describes the parabolic shape shown also in sketch (s). It 

is evident from sketch (s) that the sum of Cmg and Ch. will be a 
minimum at some value of the static margin, and that the sign 

of Cmg + Cmg at that point determines whether or not a region of axis 
positions will exist over which the wing can experience negatively 
damped rotary oscillations. These qualitative statements may be written 


explicitly by considering the equation for the damping in pitch about an 
arbitrary axis, 


Cmg + Cm = [oma | + [ema | = “ee [es | + "ns. - 


O 


2 EJ Cr, (26) 


where again Ax, refers to the distance of the center of gravity fre. 
the aerodynamic center, and the subscripted terms are referred to an 
axis through the aerodynamic center. 
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Taking the derivative of equation (26) with respect to AM , and 
setting the result equal to zero, there is obtained the axis location 
at which. the damping in pitch is a minimum 


Sm. { [era] + Qu] - W cry = 
O 


a[Axo/€] 
* [Chg] * Cra 
C 


1010. 


When equation (27) is inserted into equation (26), the minimum value of 
the damping in pitch is given as 


{epe } IE F 


ood 
min. O 


and hence, a region of instability will exist if 7 


{ [oma | + [os | 22 


If equation (29) 1s greater than zero, the boundaries of the region of 
axis positions over which instability is possible are of course given by 
setting equation (26) equal to zero and solving for AC. 


(27) 


e * Cle 0 (29) 


— Cni [Cm J, + (Ong, I. 
A 9 20 * 


Notice in both sketch (s) and equation (30) that for & given Mach number 
there will be two axis positions at which the damping in pitch vanishes. 
Then if the above procedure is carried out for a series of Mach numbers, 
one may trace out a curve as shown in sketch (t) which forms the locus 
of Mach numbers and exis positions at which the damping in pitch is zero. 
This locus thus delineates the regions of Mach number and axis position 
for which dynamic instability is and is not theoretically possible. 

Such loci, covering a wide variety of wings at supersonic speeds, have 
been presented by a number of authors. Watkins, for example, presents 
supersonic boundary curves for rectangular and triangular wings in refer- 
ences 16 and 17. At subsonic speeds, Miles! reduction of Possio's 
development to first order in frequency (ref. 11) can be used to form a 
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stability boundary curve for the two- 
dimensional wing for a given (small) 
reduced frequency. Such a curve is pre- 


sented for the entire Mach number range AX, 
in figure l, where, here + x, is the C 


distance of the axis of rotation behind 

the leading edge and k, the reduced fre- 
quency, is 0.011 for subsonic speeds and 
approaches zero for supersonic speeds. 

Notice in figure 1 that at both subsonic 

and supersonic speeds, the range of Mach 
numbers over which dynamic instability is o 
possible is largest for center-of-gravity 4.0 M, 
positions forward of the aerodynamic 
center. Further, the largest range of 
axis positions over which dynamic insta- 
bility is possible occurs near Mg = l. Sketch (t) 

Both of these characteristics are true as well for three-dimensional 
wings at supersonic speeds (see, e.g., refs. 16 to 19). 


stable 


unstable 


Effect of Mach number.- Next consider the effect of Mach number on 
the damping in pitch of a two-dimensional wing with axis at the leading 
edge. The variation with Mach number of the indicial pitching-moment 
response to a change in a will first be examined, using the information 
given in the previous sections and the indicial curves given in refer- 


ence 4. At supersonic speeds, the manner in which ~= oF (ꝙ) d, the 
2 O 3 


area corresponding to Cm; is affected can then be assessed and com- 
pared with Cmolda). At subsonic speeds, use is made of equation (18b). 
it is evident du equation (18b) that by fixing k and choosing ($4, 


such that the quantity z G(q,,k) is the same at each Mach number, one 


is free to compare finite areas - L 5 FS () dꝰ on an equivalent basis. 


As has been mentioned previously, the starting lift, at any Mach 
number, is ho,/Mo and is concentrated at the midchord. At Mo = O, 
therefore, there is an initial infinite pulse in the pitching moment 
about an axis coincident with the leading edge after which the indicial 
curve drops to /t and begins to grow asympotically toward its steady- 
state value x/2. At low subsonic Mach numbers other than zero, the 
initial pitching moment is less than infinite but very large, being 
2/Mo „ and then falls before growing toward the steady value x/2ß8. As 
the Mach number increases toward 1.0, the starting pitching moment 
falls while the asymptotic value grows, until at io = 1.0 the indicial 
curve becomes unbounded in asymptotic moment. 
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As seen in sketch (u), the effect of increasing the Mach number at 
subsonic speeds is therefore to increase rapidly the area corresponding 


to the destabilizing moment contribution, - i" F, (O) d. 
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In sketch (u), k bie chosen to be 0.011 and the values of 9, were 
picked such that - = + G(9, „k) was +4.88 at each Mach number. In the 


following discussion, the damping moment - £ Ya Fg (d + 4.88 will be 
referred to as Cms for convenience. 


At supersonic speeds, the initial value of the pitching moment 2/Mo 
continues to drop with increasing Mach number, but bere the steady-state 
pitching moment also begins to fall and at a faster rate than the start- 
ing moment, being 2/8. Even more important, as the Mach number increases, 
the number of half-chord lengths traveled to reach steady state decreases 
rapidly, being 22, for example, at Mo = l.l, as compared to 4 at Mo 2. 
As seen in sketch (v), the area representing Cm; therefore shrinks 
rapidly with increasing supersonic Mach number and becomes relatively 
unimportant at Mach numbers greater than 2. The trend of Cm with 
Mach number through the range O<M,<2 is more clearly evident in 
sketch (w). It is seen that Cms is positive, or destabilizing, 
throughout the Mach number range and that its effect is most important 
at Mach numbers near 1.0. Also shown plotted for comparison in sketch (w) 
is the variation of Cmg With Mach number. When the parameters are 
added, it is evident that the damping moment Cn, + Cms for the two- 


dimensional wing with axis at the leading edge is destabilizing in the 
Mach number range O<M,<1.414. 
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Effect of aspect 
ratio. To illustrate the 
effect of aspect ratio, 
it is convenient to 
compare the supersonic 
damping-in-pitch char- 
acteristics of a group of 
triangular wings having 
subsonic leading edges. 
The wings are of equal 
area and differ only in 
aspect ratio. As was 
done previously, the 
indicial lift responses 
to a change in angle of 
attack will first be 
examined. The effect of 
aspect ratio on the char- 
acteristic area repre- 
senting Cz. can then be Sketch (w) 
assessed. 


As has been mentioned previously, the starting lift coefficient 
after a step change in a is independent of aspect ratio and is there- 
fore equal to UG, for each wing. The parameter og, the number of 
half M.A.C. lengths required to reach steady state, is also the same for 
each wing, being a function only of Mach number. The steady-state lift 
coefficient, on the other hand, is a function of aspect ratio and 
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with aspect ratio (ref. 20). It is apparent that although C becomes 


more stabilizing with increasing aspect ratio, the destabilizing effect 
of Cms predominates, and the trend of the net damping moment is seen 
to become highly destabilizing as the aspect ratio is increased. 


By the same reasoning, the variation with aspect ratio of the 
damping moments of other types of wings can be shown to be similar (see, 
e.g., refs. 11, 21, 22, and 23). A notable exception, however, is the 
triangular wing with supersonic leading edges, whose damping in pitch 
has been shown to be independent of aspect ratio (see refs. 18 and 19). 
This characteristic may be anticipated from a study of the indicial 
response curves, since not only are the initial pitching moment Cm, (o) 
and the half M.A.C. lengths traveled to reach steady state (oa) inde- 
pendent of aspect ratio, but, unlike the subsonic-edged triangular wing, 


the steady-state p T e moment Cy (oa) is also independent of aspect 
ratio, being N „ Inspection of the results of reference 4 then 
e 


reveals that the indicial variation Cg (9) between zero and og and 
the steady-state parameter Cg (ca) are likewise independent of aspect 
ratio. 


Effect of plan-form shape.- Next, consider the effect of plan-form 
shape on the supersonic damping in pitch of a group of wings having the 
same aspect ratio. In order to make use of readily available theoretical 
results and yet compare these results with the experimental data to be 
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given later, three wings of aspect ratio 3 are chosen, having the trian- 
gular, swept, and rectangular plan forms shown in sketch (z). 


Sketch (z) 


As has been shown in the section entitled "Effect of static margin, 
the tendencies of the wings toward dynamic instability can be compared 
comprehensively by plotting their stability boundaries. For this com- 
parison, then, use is made of equation (30). The stability derivatives 
which appear in equation (30) were computed from the theoretical results 
of references 11, 20, 22, and 23. Results of. these calculations are 
shown in figure 2, wherein the stability boundaries for the three wings 
are shown as a function of axis position and Mach number. (Note that 
the axis position for each wing is measured as the distance from the 
leading edge of the M.A.C. of the wing, and that the dimensions are non- 
dimensionalized on an equivalent basis by referring them to the M.A.C. 
of the triangular wing.) It is clear from inspection of figure 2 that 
at any Mach number the triangular wing has the smallest range of axis 
positions over which dynamic instability is possible and the rectangular 
wing, the largest. 


The differences in the damping characteristics of the triangular 
and rectangular wings will be more clearly understood by a qualitative 
study of' their indicial lift responses for a Mach number of 1.2, and an 
examination of the distribution of loading due to & for the two wings. 
Consider first the indicial lift responses. 


Again, the starting lift coefficient is independent of plan-form 
shape and is l, for each wing. For the rectangular wing, the 
Lift drops abruptly after time.zero due to the loss in lift in the 
regions of the wing influenced by the formation. of the tip Mach cones 
and the starting waves from the side edges (see sketch (m)). ‘Then as 
the wing begins to emerge from the influence of the starting waves, the 
lift begins to recover, rises toward its steady-state value (given 
by eq. 6.3-2 of ref. 11), and attains this value after 12 half-chord 
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lengths of travel (eq. (24)). The variation is shown in reference 4 and 
is redrawn in sketch (aa). As mentioned previously, theoretical indicial 
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Sketch (aa) 


lift results have not yet been developed for triangular wings having 
subsonic leading edges. However, the variation shown in sketch (aa) is 
considered to be a reasonable estimate of the true shape, being based on 
knowledge of the steady-state lift (ref. 20), the time to reach steady 
state (eq. (24)), and the assumption that the shape of the variation 
would be similar to that of the wide triangular wing (ref. hà), The 

curve was adjusted within the known end points until the area correspond- 
ing to Cre agreed with that given for this parameter in reference 20. 
It is evident from examination of sketch (aa) that because of the initial 
loss in lift and the larger steady-state lift for the rectangular wing, 
Cry for this wing is significantly more negative than that for the tri- 
angular wing. Next, it is shown in references 11 and 21 that with the 
exception of regions influenced by tip Mach cones, the loading due to a 
(for Ọ >o,) for wings having swept-back leading edges increases linearly 
from zero along rays from the apex; whereas for rectangular wings the 
loading due to d increases linearly from zero along chord lines. These 


characteristics place the center of loading due to G. dr Co for the 


triangular wing and approximately z Co for the rectangular wing.” 


Then for an axis of rotation passing through the aerodynamic centers of 
the wings, the moment arm for the lift due to a for the triangular 


TDue to the influence of the tip Mach cones, the center of loading due 
to d. is shifted forward somewhat from the position it has for the 
two-dimensional wing. Calculations for the A = 3 rectangular wing 
at Mj = 1.2 show that the center of loading is at 0.605 cg. 


. 
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wing is = Co or 8 c; whereas for the rectangular wing it is approxi- 


mately z Go Thus, not only is the negative out-of-phase lift contri- 


bution CLs for the rectangular wing significantly larger Span that for 
the triangular wing, but the destabilizing damping moment - = Chg is 


larger yet, due to the larger moment arm. Calculations for the steady 
pitching parameter for an axis through the aerodynamic center ö 
(refs. 11 and 20) then reveal that C for the triangular wing is more 
negative than for the rectangular wing. The net result is therefore & 
considerably larger damping moment for the triangular wing than for the 
rectangular wing. The result of this comparison, however, should not be 
interpreted as a recommendation that the triangular rather than the rec- 
tangular wing be used on aircraft from a dynamic stability standpoint. © 
To obtain adequate static stability, the rectangular wing would generally 
be employed in combination with & tail surface, whereas the triangular 
plan form may be sufficiently airworthy without the use of a tail. The 
addition of a tail surface in effect reduces the aspect ratio of the 
rectangular wing, which reduction, as noted previously, has a highly 
stabilizing effect on the damping in pitch. The tailless triangular wing 
may therefore experience more difficulty at Mach numbers near 1.0 than a 
rectangular wing-tail combination. 


Effect of frequency.- The previous discussion has been restricted 
to the analysis.of a harmonic motion that is of vanishingly small fre- 
quency. This limitation arose as a consequence of discarding all but 
first order in frequency terms in the expansions of equations (12) 
and (14). The question arises: when the frequency can no longer be con- 
sidered small, what effect has the frequency on the damping in pitch? 


Previously, the trigonometric terms in equation (144) were expanded 
and, assuming k to be very small, terms of order k= and higher were 
eliminated. It was then found that the loss in damping from that pro- 
vided by the steady damping parameter Cna( 0a) was associated with the 
a E moment contribution corresponding to the term 


RI Fg () do. Now, however, we discard the restriction of small k 
cu perform graphically the integrations evident in equation (14) for 
supersonic speeds for several values of k. The procedure is indicated 
in sketch (bb). It is apparent from sketch (bb) that the effect of 
increasing k is to reduce the area corresponding to the destabilizing 


2 
moment contribution f x F,( 9) cos Kd. There appears another 
2 


Sagan, due to the influence of tip Mach cones, the aerodynamic center is 


shifted forward from = cg to 0.443 eg. The moment arm is therefore 
0.163 co. 


s 
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Sketch (bb) 


destabilizing contribution, k [^ F. (o) sin kPa), but quite evidently 
its effect is small compared to the reduction in the term 


-f de p (9) cos kpap. Notice further in sketch (bb) that the effect of 
increasing k becomes of marked importance when the half-period of 
oscillation is the same order of magnitude as the time for the indicial 
response to reach steady state. As shown in sketch (bb) for the fre- 
quency ks, the area ( then begins to subtract from o so that the 


destabilizing contribution - k d F,(9) cos kd can be very small. We 
may therefore expect that increasing the frequency of oscillation has a 
stabilizing effect on the damping in pitch. This conjecture is sub- 
stantiated in figures 3 and 4, where the supersonic stability boundary 
curves for aspect ratio 4 triangular and rectangular wings are shown 
plotted for various reduced frequencies. These curves were obtained 
from calculations based on the results of references 16 and 17. Notice 
that for both wings the region of possible instability is diminished as 
the frequency is increased. 


From the results of the analysis for supersonic speeds, we may 
further expect that the stabilizing effect of increasing the frequency 
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will be of even more importance at subsonic speeds, for here the indicial 
variation Fg(?) dies out at infinity. The half-period of oscillation 
is therefore always smaller than the time to reach steady state. 


The situation for the term - f^ Fa(9) cos k9d? is shown in sketch (cc). 
O 
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Sketch (cc) 


It is evident in sketch (cc) that the destabilizing moment 
- L Fg (O) cos kPa? diminishes rapidly as the frequency is increased. 


The effect of this reduction on the damping in pitch can be illustrated | 
by plotting the subsonic damping-moment coefficient against reduced 
frequency for the two-dimensional wing (with axis at the leading edge) 
for Mach numbers O, 0.5, and 0.8. The results, which were obtained 

from reference 4, are shown in figure 5. It is seen in figure 5 that 
the garge destabilizing effects of the moment contributions 


— f. " Fa(9) cos Kk and K S F (C) sin k~aP are confined to a rela- 
tively narrow range of reduced frequencies. Notice further in figure 5 
that the range of frequencies for which instability is possible is small 
at Mo = O (0<k<0.04) and grows with increasing Mach number. This is 
believed to be the primary reason why unsteady lift effects were found 
to be unimportant at low speeds but are of great importance at speeds 
near the sonic speed. 


. 
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Effect of thickness.- In this section, the concepts regarding the 
indicial functions discussed in previous paragraphs will be used to 
examine the effect of profile thickness on the damping in pitch of & 
two-dimensional wing flying at a supersonic Mach number. 


It will be recalled that for low frequencies the damping in pitch 
may be considered to be the sum of the damping moment due to steady 
pitching, Ong oa), and the area representing Cm; on the plot of the 
pitching-moment response to a step change in angle of attack 
(sketch (k-ii)). The airfoil thickness, of course, influences both 
these parameters. For the purpose of this analysis, it will be con- 
sidered sufficiently accurate to study the effects of thickness only on 
the steady parameter (og) and on the end values Cy,(0), Cn (oa) 
and dg of the indicial response curve. It will then be assumed that 
shape of the indicial response curve is not appreciably altered by thick- 
ness effects, whereupon the major effect of thickness on Cms can be 
assessed by adjusting the response curve given by the linearized theory 
to fit the corrected values of Cy (0), Cm(og) and og. 


It has been shown by Busemann (see, e.g., ref. 2h) that the loading 
at an element of a lifting airfoil may be represented by 


= = Ci [au- za. - Cs [orit [0407] } (31) 


where 
2 


J Mor -1 
 [ru]Wo* - h[w5*-1] 
| 21M 2-115 


C1 * 


02 


and qu and c, are the local slopes of the upper and lower surfaces, 
respectively, measured with respect to the chord line, For the special 
case of an airfoil that is symmetrical about the chord line, 

c, equals c} and equation (31) reduces to the following: 


OP . 200, + 10020 (32) 
do 1 Eu 


For the sake of simplicity, in what follows equation (32) will be used 
for the local loading rather than equation (31). The results will there- 
fore be valid only for airfoils that sre symmetrical about the chord line. 


———É 


42 ONF IDENTIAD NACA RM A52LOla 


The lift on the airfoil due to angle of attack is obtained simply by 
integration of equation (32) across the chord. 


B OD 
CL (oa) Co Í do dx 
C hos po 
CLa (oa) = aie = 2C4 + = L 7 o, dx ( 33) 


Similarly, the moment due to angle of attack, measured about the leading 
edge, is | 


—Ó [° & 
Cn We pe 
CN ee oe o - [o eux ax (30) 


When both the leading and trailing edges are on the chord line, equa- 
tions (33) and (34) reduce to the following: 


2C1 


Cra 08) 
(35) 

2RCe 

2 


cmd, (ca) = = C1 + 


Co 
where Q is the enclosed area of the airfoil section. 


The local loading given by equation (32) is still applicable for 
the steady pitching case with this added provision: that now a repre- 
sents the local angle of attack. For an airfoil pitching about its 
leading edge, the variation of a with chordwise distance is given 
by a = qx/Vg, where q is the pitching velocity. Then again, for the 
Lift due to pitching, 


|r 

m 

eB 
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Cr, (oa) = 


| 
, 
Jk 


Cha! (c4) = qu x dx (36) 
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Likewise, the pitching moment due to pitching about the leading edge is 


1 Co 
Cm () - — = 
In a) co? 55 


Cm' (04) 4 8Cə co 
„ a Data 


The primes on the parameters indicate that the moments are measured about 
and the wing is pitching &bout the leading edge. 


In order to illustrate the result given by equation (37), the 
damping moment due to steady pitching has been computed for a family of 
airfoils having symmetrical biconvex parabolic arc sections. The 
results are presented in figure 6. The dashed line in figure 6 repre- 
sents the theoretical locus of Mach number and leading-edge angle for 
bow-wave detachment (ref. 24). Since the requirements of the second- 
order theory are violated to the left of this line, the damping curves 
there are not rigorously correct, and therefore are terminated a short 
distance beyond the line. 


To determine the effect of thickness on the parameter, Cms', it is 
assumed that the shape of the indicial pitching-moment variation is not 
significantly affected by thickness. When this assumption is made, it 
becomes necessary only to correct the initial and final ordinates of the 
indicial curve Ci (o) and Cm (o. and the number of half-chord lengths 
traveled to reach steady state, og. The response curve is adjusted to 
fit these new end values, so that the area representing the cor- 
rected Cys' can then be determined. 


The effect of thickness on the final ordinate of the indicial curve, 
'(cg), has already been determined (eq. (34)). Now consider the 
Ead lift and moment. 


Assume first a thin flat wing to be flying at zero angle of attack 
in a uniform stream of density jpg. At time zero, it starts to sink 
with downward velocity voc. The impulsive start causes a plane com- 
pression wave to be emitted from the lower surface at the speed of sound, 
which is constant throughout the stream. (The pressure and therefore 
the temperature is sensibly constant.) At the same instant, a plane 
expansion wave is emitted from the upper surface and travels at the same 
velocity in the opposite direction. At the end of an infinitesimal time 
period At, the wing has sunk a distance Voc. The sound waves have 
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moved a distance a ,At in both directions. The mass of fluid affected 
is 2pgeg&8ogAt. Then, from Newton's second law, one may write 


(fy-f2) At = 2pgegag^t Voa 


where f, and f; are the forces per unit span on the upper and lower 
surfaces, respectively. Converting to coefficient form, it then appears 
that 


For the wing of finite thickness for which the pressure varies along the 
chord, a different result is expected. The wing is assumed to have been 
flying at zero angle of attack for a period of time long enough for the 
flow to have established itself around the surface. In this case, the 
density, temperature, and therefore the speed of sound, vary along the 
chord due to the curvature of the surface. At time zero the wing again 
starts to sink, but the compression and expansion waves are propagated 
at different rates at various chordwise positions so that 


At [fy-f1) (x) = 2 Um (x) E. (x) ot Voa 


and the starting lift becomes 


cr (o) h Co p & 4. | 
= —— — dx 8 
Q MoCo f Po &o at 


From reference 24, the variation of o/o a/ag with Cp is found to be 


„ 2 
D a 7 27 
c = E ST e, | (39) 


where Cp is the pressure coefficient at a point on the surface when 
the airfoil is at zero angle of attack. Again, from reference 2h, 
P-7Po 


Cp = 


= C40y + ese (40) 


where p and pg are, respectively, the pressure at a point on the sur- 
face of the airfoil and the free-stream static pressure. 
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For the symmetrical biconvex parabolic are airfoil considered 
earlier, the local slope qu is given by 


cy = 25 [1 - | (41) 


where 8 is the thickness ratio. 
Then expanding equation (39) in a binomial series, retaining terms 


through the second order in oy, and performing the integration indicated 
by equation (38), there results for the starting lift coefficient 


Cna(o) = 5 140. t," Lecce (42) 


A similar procedure can be carried out for the starting pitching-moment 
coefficient, measured about the leading edge 


Cm," (o) = = 


Ch Co) * 140.0hM [200,87 M 810018 ] (43) 


lt is necessary also to consider the effect of thickness on the time 
required for the loading on the airfoil subjected to a step change in 
angle of attack to reach steady state. As was seen in sketches (n) 

and (p), the time to reach steady state is exactly equivalent to the 
time required for the wing to escape the influence of the sound waves 
emitted from the leading edge at the start of the motion. For purposes 
of analysis, the situation is reversed by requiring that the wing remain 
stationary in a fluid moving with supersonic velocity Vo, and causing 
& sound wave to be emitted at the leading edge of the wing. To an 
observer standing on the upper surface of the airfoil, there will appear 
to be two wave fronts whose velocities tangent to the surface are 

[V+a] and [V-a], where V anda are the local stream velocity and local 
speed of sound, respectively. Then the time required for the slower 
wave front to clear the airfoil is given by 


t Co 2 
wf mf Rm Len. (uz) 
5 o  [V-a]cos[tan "o4] Vo o IV/ao-/ ao! 
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It is therefore necessary to obtain V/ag and a/ao as functions of 
chordwise location x. Again, in reference 2h these quantities are 
found to be functions of the local pressure coefficient Cp · 


y-1 
2 
Eu E 4 rMo op | e 
80 2 
(45) 
5 4. 
S us xit is peme) 
Bo 74 80 . 
Expanding equations (45) through the second power of Cp and Vu 
through the second power of o,, there is obtained 
" co l + 8 ou? | dx 
m (o ce chc m (46) 
Vg [M5-1.] - l-k CPHR Cp 
where 
K. = 80 + O.10Mo* 
1 Mo - 1 
"e 0.125M5[Mo9?-1] + 0.03M,* 
2 Mo d 1 


The pressure coefficient Cp its now assumed to be that existing at 
a point on the surface when the airfoil is at zero angle of attack. This 
assumption, of course, introduces some error since the angle of attack 
has changed discontinuously from zero to a at t= O. However, if we 
require that a be less than 8*, then the error introduced through 
neglect of the effect of the change in a on the pressure coefficient 
is negligible. Then, as in the starting lift problem, 


Cy = C404 + Cao? (4-7) 


Inserting equation (47). into equation (46), dividing numerator by denom- 
inator, and again retaining terms through the second power in cu, there 
results l 
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For the symmetrical biconvex parabolic are airfoil, Oo, is replaced by 
equation (41), whereupon for this airfoil the time to reach steady 
state becomes 


ta = oto [14£] (49) 
Vo[Mo-1] 


where 


Oa, = Zo ta "C 11751 (20). 


The effect of thickness on the damping parameter Cm! can now be 
estimated approximately by &djusting the indicial pitching-moment curve 
obtained from the linearized theory to the end values given by equa- 
tions (34), (43), and (50), and then measuring the area corresponding to 
the "corrected" Cms- The process is indicated in sketch (dd) below, 


and the adjustment formulas are given by equations (51). 


. 
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Sketch (dd) 
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Cm. (9, ) - Cn (Po) = md (0) = — |Ana" (oa) = Ang! (9) 
i : (51) i 


where 
Ana’ (0) On, (e) = Cmo, ' (0) 
mer (7a) = Cmq. (Co z nd. (dag) 


and the subscripts o and i indicate values corresponding to the linear- 
ized theory and the adjusted theory, respectively. Notice also in 
sketch (dd) that the area corresponding to the adjusted value of Cms. 
is given by the following expression: 


Cms! = Cm ' (O- oa, Any! (ca) (52) 


. The parameter Cms' for pitching about the leading edge has been 
computed, using equations (35) to (52) for the same family of airfoils 
as that discussed previously and the results added to those obtained 
for Cg (ea. (37)). The total damping moment Cm,' + Cm, is shown 
plotted as a function of Mach number in figure T. “te is evident from 
examination of figure 7 that the effects of thickness are small and are 
in the destabilizing direction with increasing thickness ratio. It 
should be noted that this result does not agree with those of Jones and 
of Wylly in references 25 and 26 which indicate that the effect of thick- 
ness is extremely stabilizing. In reference 27, however, Van Dyke takes 
issue with the results of references 25 and 26, and offers an alternate 
solution correct to the second order in thickness. From this result, 
the following relationship for the damping-moment coefficient for pitch- 
ing about the leading edge can be deduced: 


ese] pe] [+ CEH] c 


where bend Cni is the damping-moment coefficient obtained from the 


linearized theory and 
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y = Ztl Mo? 
2 p? 
Ace for biconvex airfoils 
À -i for double-wedge airfoils 


The trend of the damping-moment coefficient with thickness obtained from 
this result agrees closely with that obtained in the present paper. 


EXPERIMENT 


Wind Tunnel 


The experimental investigation of the damping-in-pitch character- 
istics of the various model configurations was conducted in the 
Ames 6- by 6-foot supersonic wind tunnel. In this tunnel, the Mach num- 
ber can be varied continuously within the subsonic speed limits 0.60 to 
choking, and at supersonic speeds from 1.20 to 2.00. The total pressure 
can also be changed at any time within the limits 3 pounds per square 
inch absolute to atmospheric. A complete description of the wind tunnel 
is given in reference 28. | 


Models 


The pertinent dimensions of the five wing-body configurations 
investigated in the present report are shown in figure 8. All the wings 
had sections which were symmetrical in streamwise planes and 3 percent 
thick. Wings having leading-edge sweep angles less than or equal to 45° 
(wings c, d, and e in fig. 8) had biconvex circular-arc sections in 
streamwise planes with maximum thickness at 50-percent chord. Wings 
having leading-edge sweep angles greater than 45° (wings a and b in 
fig. 8) had NACA 0003-63 sections in streamwise planes. The wings were 
identical in plan form and section shape to those investigated in the 
series of force tests reported in references 29 to 33. Strength consid- 
erations required that wings c, d, and e be built of steel. Wings a 
and b were built of aluminum. 


The models were fitted with bodies such that the distance from the 
apex of the body to the wing-body intersection was the same for each 
model. Also, the ratio of wing area to body maximum cross-sectional area 
was 17.9 for each model. Forward of the point of maximum radius, the 
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bodies were of thin laminated wood construction and were identical in * 
shape to those of references 29 to 33. The aluminum afterbodies were TE 
cylindrical in shape and terminated at the trailing edges of the wings 
to permit the models to be deflected with sufficient amplitude when 
mounted on the support structure. 


Model Support System 


The support system used in the present investigation was the same : 
as that developed for the single-degree-of-freedom free-oscillation 
experiments of reference 10. The interested.reader is referred to that 
report for a detailed description of the oscillation mechanism. For —- 
the present tests, two changes were mađe in the system which served to 
improve the accuracy of the data. The changes and the reasons for them 
are as follows: 


l. The vertical flexure pivots which provided vertical restraint 
in the tests of reference 10 were replaced by bearings after it was 
found that the heavier models of this investigation caused the pivots 
to twist laterally when the model was given its initial displacement. 
An undesirable yawing oscillation was thereby induced which destroyed 
the linearity of the system. Installation of the bearings entirely elim- 
inated this yawing tendency. The damping due to friction in the bearings s 
was somewhat larger than that of the flexure pivots, but by frequent = 
checking of the wind-off damping and regular replacement of the bearings 
it was possible to maintain the tare damping well within acceptable - 
limits. 


2. A stiffening strut was installed between the sting support and 
the tunnel ceiling in order to remove the possibility of coupling between 
the model oscillation and the resonant mode of the sting support. The 
influence of the strut on the aerodynamic forces at subsonic speeds was 
investigated by removing the strut and recording data for one of the 
models at a frequency sufficiently below the sting resonant frequency to 
avoid excessive sting vibrations. Results obtained with the strut 
installed and removed agreed within the experimental accuracy. 


By effectively eliminating both yawing tendencies of the model and 
vibrations of the sting support, it can be seen from comparison of the 
present results with those of reference 10 that the deviations between a 
number of observations at a given Mach number were markedly reduced. 
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Scope of Tests 


Investigation of the damping-in-pitch characteristics of the various 
model configurations was conducted over a supersonic Mach number range 
of 1.20 to 1.90 and, where possible, over a subsonic Mach number range 
of 0.60 to 0.90. The Reynolds number for the major portion of the tests 
was held constant at 1.6 x 108 per foot. Additional data for one of the 
models was obtained at Reynolds numbers of 1.0 x 108 per foot and 
3.2 x 108 per foot. 


All the tests were conducted with the models &t & mesn &ngle of 
attack of 0°, the angle being measured from the mean line of the sting 
support to the axis of the test section. For each oscillation, an ini- 
tial displacement of +5° was imparted to the model by means of the 
pneumatically actuated pawl arrangement described in reference 10. 


Other important variables are listed below for the five wing-body 
configurations investigated: 


Range of moment Axes of 
of inertia, rotation, Range of reduced 
Model slug - ft^ $ M.A.C. frequency 

A = 2 triangular 0.0534 - 0.0589 35, 45 0.012 - 0.037 
A = 3 triangular .0368 ~ .0ho2 25, 35 .018 - .050 
A = 4 triangular 0h91 - .0512 25, 35, 45 .0O0T - .ohi 
A = 3 swept -0549 - .0576 25, 35 OIL. .025 
A = 3 unswept .0h22 - .0h88 20, 35, ho .OO0T - .025 


Reduction of Dat& 


The technique used in this investigation for reducing an experi- 
imental oscillation-decay record to the damping-in-pitch coefficient 
Cm, + Cm, was identical to that described in reference 10. Briefly, 
this technique involves plotting the envelope of the oscillation-decay 
curve against time on semilogarithmic graph paper. If the damping of 
the system is & linear function of the angular velocity, this plot will 
be a straight line, the slope of which is proportional to the damping 
term. Calibrations for the model moment of inertia and the damping due 
to internal friction, and standard measurements for the density and 
velocity of the air stream then enable one to reduce the aerodynamic 
damping term to coefficient form. The aerodynamic restoring moment Cm. 
can be obtained from measurements of the frequency of oscillation and 
& calibration for the static spring constant of the system. Results for 
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the restoring moment about two axes of rotation a known distance apart 
then permit calculation of the lift-curve slope CL: 


Corrections to Data 


No wind-tunnel-wall corrections were made to the subsonic results 
for any of the measured stability derivatives. To the author's knowl- 
edge, no corrections have been developed for application to the forces 
acting on a finite-span wing oscillating in & compressible fluid. 
Tunnel-resonant-frequency effects were investigated by use of refer- 
ence 34, and it was found that the range of frequencies used in this 
investigation was well below the lowest calculated tunnel resonant fre- 
quency. In view of the fact that the subsonic damping moments measured 
for the A=4 triangular wing in the present investigation agreed 
reasonably well with the results for a similar wing obtained during an 
investigation in the Ames 12-foot wind tunnel (ref. 35), it is believed 
that the effect of the tunnel walls on the SAID in piteh was not 
significant. 


A subsonic correction could have been made to the static parameters 
Cig and Cy,- However, since this correction would have been very small, 


and in view of the uncertainty involved in applying static corrections 
to the results of EE measurements, it was decided not to make the 
. correction. 


The effects of constriction of the flow &t subsonic speeds by the 
tunnel walls were taken into account by the method of reference 36. 
This correction amounted to, at most, a 2-percent increase in the Mach 
number &nd in the dynamic pressure over that determined from & calibra- 
tion of the wind tunnel without a model in place. 


For the tests at supersonic speeds, no corrections were required to 
be made on either the aerodynamic measurements or the air-stream 
properties. 


The reader is referred to reference 10 for a discussion of the 
precision of the data. Since the method of reduction of the data was 
the same as that of reference 10, the uncertainty in the recorded value 
of a given damping coefficient remained of the order of £0.02. However, 
by virtue of the elimination of yawing tendencies in the model and sting- 
support vibrations, the standard deviation of a number of observations 
at a given Mach number was reduced from £0.06 to 40. 03. 
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RESULTS 


Damping Coefficients 


The basic experimental results for the damping coefficient Cmq + ing, 
for the five wing-body configurations investigated are shown in fig- 
ures 9 to 13. In cases where data could be obtained at subsonic as well 
@s supersonic speeds, these results are shown plotted through the Mach 
number range 0.60 to 1.90 for ready comparison of the magnitudes of the 
coefficients in the two speed ranges. Also shown plotted in figures 9 
to 13 are the theoretical results for the damping coefficient Cm, + Cm. 
at supersonic speeds for the wings alone. In addition, theoretic 
results for the triangular wing-body combinations, considering inter- 
ference effects, are shown in figures 9 to ll. The sources of these 
theoretical results are references 11, 16 through 23, and 37. The exper- 
imental results for subsonic speeds are not compared with theory, since 
& rigorous linearized theory for triangular, swept, or unswept wings 
oscillating in a subsonic compressible flow has not yet been developed. 
A very approximate analysis for Cm, and Cm, for the aspect ratio 4 
triangular wing was presented in reference lO. In this development, the 
parameter Cg. was approximated using the incompressible theory of ref- 


erence 8 for the elliptic wing. However, a more rigorous analysis for 
the elliptic wing in compressible flow (ref. 11) has indicated that the 
result of reference 10 for Cm may be incorrect. The good agreement 
between the experimental subsonic results of reference 10 and the sub- 
sonic theory of reference 10 may be fortuitous. 


The results for the three triangular wings (figs. 9 to 11) and the 
swept wing (fig. 12) show that in the supersonic speed range, in every 
case, the Linearized first-order-in-frequency theory provides a reliable 
guide for obtaining both the magnitude and trend with Mach number of the 
damping coefficients. Of particular importance are the theoretical 
predictions of ranges of Mach number over which the wings could experi- 
ence negatively damped oscillations. The experimental results for 
the A= 3 triangular wing (fig. 10) are inconclusive in this regard 
Since the theoretically predicted unstable range occurs in the low super- 
sonic Mach number range, 1.0 to about 1.08, where it was not possible to 
obtain data. The predictions by the theory of larger unstable ranges 
for the A s h triangular wing and the A= 3 swept wing were borne out 
by the experimental results, as shown in figures 11 and 12 by the exper- 
imental points plotted below the abcissas. In the absence of available 
theoretical results for unswept wings, the experimental results for the 
unswept wing are compared with theoretical damping coefficients for a 
rectangular wing of aspect ratio 3 (fig. 13). It is seen from examina- 
tion of figure 13 that the theoretical prediction of a rather large 
region of instability for an axis of rotation at O. 208 (fig. 13(a)) is 
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borne out by the experimental results but that the theory fails to pre- 
dict the reversal in trend of the damping coefficients at low supersonic 
Mach numbers for axes of rotation closer to the aerodynamic center 
(figs. 13(b) and (c)). It is clear from these results that the effects 
of taper can be important and should be considered for an accurate theo- 
retical appraisal of the damping-in-pitch characteristics of unswept 
wings. 


It should be mentioned that at a Mach number of 0.90, the damping 
of the A = triangular wing (fig. 11) and the A = 3 triangular and 
swept wings (figs. 10 and 12) were highly nonlinear with angle of attack, 
being stabilizing for angles of attack of 5O and decreasing steadily to 
zero at about 1°. The dashed lines for this Mach number shown in fig- 
ures 10, 11, and 12 are intended to indicate the range of this nonlinear 
variation. Further discussion of the phenomenon is withheld to a later 
section of this report. The above-mentioned difficulty was not encoun- 
tered with the A = 2 triangular wing, and the value of Cmq + Cmá 
at Mo = 0.9 for this wing (fig. 9) represents the damping throughout 
the amplitude range. Subsonic results for the unswept wing could not be 
obtained due to the highly erratic behavior of the aerodynamic restoring 
moment which generally was so strongly destabilizing as to counterbalance 
the spring restoring moment and force the model against its stops. 


Reynolds Number Effects 


For the A = 2 triengular wing at supersonic speeds (fig. 9), 
damping coefficients have been obtained at three Reynolds numbers rang- 
ing from 1.18 x 108 to 3.77 x 10°, while at subsonic speeds (fig. 9) they 
have been obtained at 1.18 x 10 and 1.89 x 108. In both cases, it is 
evident that at least within the range of Reynolds numbers at which data 
were recorded there is no significant effect of scale. It is reasonable 
to assume that a similar lack of scale effects exists for the other wings 
ag well. 


Aeroelaetic Effects 


In figures 14 to 23, the static parameters Cm, and CL, for the 
five wings are compared with force-test measurements of similar wings 
(refs. 29 to 33). Such & comparison is useful as an indication of the 
importance of aeroelastic effects on the aerodynamic properties of the 
less rigid models of the present tests. The fact that at supersonic . 
speeds there is excellent agreement between the two experiments for the 
lift-curve slope CL? whereas the present results for Cm, are 
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consistently smaller,is attributed primarily to the difference in the 
body shapes of the models of the two investigations. In the force-test 
investigation, the afterbodies were carried beyond the trailing edges of 
the wings, which would tend to shift the centers of pressure of the wing- 
body combinations farther rearward while leaving the total lift essen- 
tially unchanged. From this reasoning, it is concluded that aeroelastic 
effects on the aerodynamic coefficients of the present tests are of 
secondary importance. At subsonic speeds, there are differences in mag- 
nitude of both CI, and Cm,- However, in view of the fact that no large 
effects of aeroelasticity were found at supersonic speeds where the 
dynamic pressure is greatest, it is believed that the differences 

in Cy, and Cy, at subsonic speeds are not due to aeroelastic effects. 
It must therefore be stated that at the present time the differences at 
subsonic speeds are not understood. 


Transfer of Axes 


Sufficient information has been given in figures 9 to 23 to permit 
the transfer of the damping-in-pitch results at supersonic speeds to 
axes of rotation other than those used during the tests. Such a trans- 
fer, if valid, of course greatly increases the range of applicability 
of the given data. The validity of the transfer equation (given in 
appendix A) depends primarily on the invariance of the stability deriv- 
atives with changes in angle of attack and angular velocity and could be 
best checked by experimental means. Checks on the method are provided 
in figures ll(c) and 13(c) for the A= 4 triangular wing and A = 3 
unswept wing, respectively. Shown in these figures are the experimental 
damping coefficients obtained during the investigation and coefficients 
for the same axis position computed, by means of the transfer equation, 
from the results of experiments at two other axis positions. It is 
evident that the transferred results portray the actual data with suffi- 
cient accuracy to establish the validity of the method. 


DISCUSSION 


Effect of Aspect Ratio 


It will be remembered from the discussion in the theoretical section 
of this report that the linearized theory indicates that a decrease in 
aspect ratio has a highly stabilizing effect on the damping in pitch of 
a triangular wing with subsonic leading edges. This indication is con- 
firmed in figure 24, wherein the damping in pitch is shown as a function 
of Mach number for the three triangular wings having their axes at 0.35c. 
It is noteworthy that the A = 2 wing shows no tendency toward dynamic 
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instability in the Mach number range 1.2 to 1.9, that for Mo- 1. 2 the 
trend of the damping coefficients of the A = 3 wing is toward insta- 
bility, and that a range of instability exists for the A - A wing. It 
is of further interest to note that for Mach numbers grester than 1.67, 
in which range the leading edges of the A = 3 and A= H wings are 
supersonic, the damping coefficients for these wings are essentially the 
same. This serves as a partial confirmation of the theoretical predic- 
tion of reference 19, wherein it is shown that the damping coefficients 


of triangular wings having supersonic leading edges are independent of 
aspect ratio. 


A more comprehensive way of comparing the tendencies of the wings 
toward dynamic instability would be to show the entire ranges of axis 
positions and Mach numbers over which dynamic instability is possible, 
that is, to plot the stability boundary curves for the three wings. 
in the Mach number range 1.2 to 1.9, however, only the A = À wing has 
a region of instability and this curve is shown, compared with the theo- 
retical curve, in figure 25. 


The preceding comparison (fig. 24) was useful primarily for verify- 
ing the theoretical prediction regarding the role played by aspect ratio 
in determining the tendency of a triangular wing toward dynamic instabil- 
ity. A comparison of the actual magnitudes of the damping moments was 
masked, however, by the fact that the coefficients for the three wings 
were referred to their own respective characteristic lengths. In the 
following comparison, an attempt is made to overcome this difficulty by 
posing the question: Given three triangular wings of aspect ratios 2, 3, 
and 4, of equal area, and required to have the same restoring moment 
(in ft-lb) at a Mach number of 1.4, how do the physical magnitudes of 
the damping moments compare? For this purpose, the aspect ratio 3 wing 
is chosen as a standard, required to have a static margin of 0.05 
at Mo =1.4, and the damping moments of all the wings referred to the 
M.A.C. of the aspect ratio 3 wing. The comparison on this basis is 
shown in figure 26. It is evident from examination of figure 26 that on 
this basis also, the effect of decreasing the aspect ratio is beneficial 
to the damping in pitch, the aspect ratio 2 wing having the largest 
damping moments throughout the range of test Mach numbers, and the Sepsey 
ratio 4 wing the smallest. 


Effect of Plan-Form Shape 


For this comparison, use is made of the results for the aspect 
ratio 3 triangular, swept, and unswept wings. First, the prediction 
made in the section entitled "Theory" regarding the tendencies of 
three similar wings toward dynamic instability is examined by comparing 
their stability boundaries. From the preceding section, it will be 


NACA RM A5PLOha — 51 


recalled, however, that no region of instability exists for the A= 3 
triangular wing in the Mach number range 1.2 to 1.9. It is apparent, 
therefore, that as predicted by the theory the triangular wing is 
superior in this regard. The stability boundary curves for the swept 

and unswept wings, with axis position referred to the leading edge of 

the M.A.C., are compared in figure 27. Again, the theoretical predic- 
tion is borne out for, as seen in figure 27, a considerably larger region 
of instability exists for the unswept wing. 


Next, the magnitudes of the damping moments &re compared by choosing 
the A 3 triangular wing as a standard, with a static margin of 0.05 
at Mo = 1.4, and, as in the previous section, requiring that the wings 
have equal restoring moments at this Mach number. The comparison is 
shown in figure 28. On this basis, it is again apparent that the trian- 
gular wing is the superior wing, followed in order by the swept and 
unswept wings. 


Notice in figures 26 and 28 that since the damping coefficients of 
all five wings investigated are referred to the M.A.C. of the A= 3 
triangular wing and since all wings have equal restoring moments at 
Mo = 1.4, the results of figures 26 and 28 may be compared directly. 
Here it appears that the A= 2 triangular wing is the most desirable 
wing from a longitudinal-dynamic-stability standpoint, and the A= 3 
unswept wing the least desirable. As mentioned previously, however, 
it should not be inferred that the result of this comparison implies 
& recommendation for the use of triangular rather than unswept wings. 
Obviously, the addition of tail surfaces could alter considerably the 
relative damping-in-pitch merits of the wings investigated. 


Effect of Thickness 


It has been shown by the results of experiments (figs. 9 to 12) 
that the supersonic linearized potential theory provides a reliable 
basis for predicting the low-frequency damping-in-pitch characteristics 
of thin wings of finite span. Considerable doubt has been shed, however, 
over the ability of the linearized theory to predict the dynamic behavior 
of wings that cannot be classified as thin. In recent years, a number 
of reports have been issued (refs. 25, 26, and 38) which indicate that 
second-order thickness effects have a profoundly stabilizing effect on 
the damping in pitch of two-dimensional wings. In particular, refer- 
ence 26 indicates that by increasing the thickness of an infinite-span 
wing from zero to only 1-1/2 percent, the rather large region of insta- 
bility predicted by the linearized theory is completely eliminated. In 
contrast to these results, the approximate analysis of this paper 
(see Theory) indicates that thickness effects are relatively small, the 


destabilizing effect on m AI canceled by & stabilizing 
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effect on Cm5. In view of these conflicting results, it is of some 
interest to obtain additional experimental data which would serve to 
clarify the issue. Of particular interest are the effects of thickness 
for finite-span wings. In this case, a limited comparison can be made 
between the experimental results for the 3-percent A = 4 triangular 
wing of this report and the results of reference 10 for two 6-percent 
triangular wings having the same aspect ratio. Aside from the differ- 
ences in profile and small differences in the dimensions of the bodies, 
the models and test conditions of the two experiments were almost iden- 
tical. The comparison is shown in figure 29. It is seen that for this 
wing, the effects of thickness are small and are in the direction of 
decreasing stability with increasing thickness. Of primary importance 
is the fact that all the wings exhibit regions of Mach number over which 
dynamic instability is experienced. 


Although admittedly the experimental evidence in figure 29 is not 
conclusive, nevertheless in the light of these results it is difficult 
to accept those results which indicate that large stabilizing effects of 
thickness &re present on infinite-span wings, especially in view of the 
fact that the behavior of a triangular wing with supersonic edges is in 
many ways quite similar to that of an infinite-span wing. 


Effect of Nonlinearities 


From examination of the experimental oscillation-decay records at 
both low subsonic and all supersonic speeds, it was -evident that the 
assumption of a linear second-order system was well justified. Aside 
from a moderate amount of scatter, a decay curve could be fitted by an 
exponential curve over a wide range of amplitudes, verifying that the 
damping parameter was essentially independent of amplitude. At high sub- 
sonic speeds, however, this linear behavior was no longer true for most 
of the wings tested. As an example, figure 30 shows the oscillation— 
decay record for the A = 3 triangular wing at Mo = 0.9%. It is seen 
that after an initial displacement the amplitude of oscillation quickly 
dies out to a low level, but there sustains itself indefinitely. 
Similar phenomena have also been reported in the forced oscillation 
experiments of reference 35 for & wing-body combination having a trian- 
gular wing of aspect ratio 4 and in flight tests of various high-speed 
aircraft, particularly tailless aircraft. The deleterious effects of 
such sustained oscillations on the qualities of an aircraft as a gun 
platform or from & structural standpoint are obvious. It is therefore 
of considerable importance to gain an understanding of the phenomenon, 
with a view toward advancing means of eliminating it. 


At the present time, insufficient information, either theoretical 
or experimental, is available to enable a complete and authoritative 
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description of the physical event; & tentative hypothesis is therefore 

advanced instead, which would appesr to fit the observed facts. It must 
be emphasized, however, that in the light of additional information, the 
assumptions and conclusions drawn here may require considerable revision. 


Specifically, the observed facts are these: (1) The damping 
coefficient is highly nonlinear, being negative or stabilizing for angles 
of attack greater than about 19, and zero at 19. The possibility that 
the nonlinearity is caused entirely by scale effects and/or wind-tunnel- 
wall interference is ruled out by the fact that similar behavior is 
observed in full-scale flight tests. (2) The frequency of oscillation 
is essentially constant throughout the amplitude range, implying that 
the restoring moment is a linear or nearly linear function of the angle 
of attack. 


From these observations, it is hypothesized that the characteristic 
differential equation governing the motion is of the form 


I& T O (,d) d + Kea = 0 


where Q(a,a&) represents the damping coefficient as a nonlinear function 
of the angle of attack and angular velocity, and I and K are constants. 
It is therefore necessary to search for a mechanism that can affect the 
damping moment to & much larger degree than it does the restoring moment. 
In this regard, it is pertinent to first review some of the character- 
istics of wings in steady transonic flow. It is well known that at high 
subsonic Mach numbers, the essential feature of the flow is the appear- 
ance of shock waves on the wing as the speed of air over the wing surface 
exceeds the local speed of sound. These shock waves cause large changes 
in the pressure and also promote flow separation near the shock wave, 

aque to the presence of adverse pressure gradients. The effects of these 
disturbances on the pitching-moment characteristics of an unswept wing 
of aspect ratio 3 are shown in reference 39 and reproduced as curve (a) 
in sketch (ee). Notice that 65 


the slope of the pitching 
moment is positive for small 
angles of attack, indicating 
that tbe &ircraft would be 
statically unstable in that 
range. It was also illus- 
trated in reference 39 that 
by reducing the curvature of 
the wing profile near the 
trailing edge, thereby | Sketch (ee) 

reducing the tendency toward separation, the positive slope through zero 
could be eliminated. This result is shown as curve (b) in sketch (ee). 
By virtue of the beneficial effect of reducing separation, it is argued 


(a) basic wing —LÓ 
(b) modified wing ———— 
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here that in the absence of separation and strong shocks on the basic 
wing, the variation of pitching moment with angle of attack would be 
nearly linear, as shown by curve (a) in sketch (ff), and that the effects 
of these disturbances can 
be lumped into an additional 
moment variation, curve (b), 
such that the sum of (a) 

and (b) then gives the 
observed nonlinear variation. 


For triangular wings, 
the reversal in slope of the 
pitching moment at small 

angles of attack does not 
Sketch (ff) appear (see refs. 29, 30, 
and 31). This, however, does not preclude the presence of the nonlinear 
variation due to separation. For a triangular wing with center of grav- 
ity at, say, the quarter M.A.C. point, the slope of the pitching-moment 
curve can be quite large. As seen in sketch (gg), the linear varia- 
tion (a) can then greatly outweigh the effect of the nonlinear curve (b), 
with the result that a scarcely evident nonlinearity appears in the 
combined curve. | 


The preceding examples were taken from the 
results of static measurements, where the angle 
of attack is simply the angle between the chord 
line of the wing and the free-stream direction. 
When the wing is oscillating, however, each 
point 5 measured from the axis of rotation 
experiences an additional angle of attack 
a&/Vo, due to the angular velocity. Then if 
the steady-state pitching-moment variation 
caused by the separation (curve (b) in 
sketch (gg)) is approximated in the range 


O<a<anax by 


Sketch (gg) 


in the unsteady case it is, 


om (a+ #) = alas] -v fas . (54) 
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where now € is assumed to be the distance from the axis to the point 
at which the additional lift due to separation is concentrated.9 Now, 
since for slow frequencies G&£/V) is much smaller than a, all but 
first-order terms in d. are neglected so that equation (54) becomes 


Cm = . [a-bo7] + = [a-3ba?] (55) 
O 


The characteristic equation of motion then becomes 


10 2 [tenga i x Le- zü [ente] \ - 0 (56) 


ev. 


where [Cmgt+Cmg] and Cm, are the (constant) stability derivatives which 
would be present in the absence of separation. Thus, it appears that 
although the effect of the nonlinearity on the restoring moment can be 
insignificant, the possibility still remains that the damping moment can 
be profoundly affected in the event that [Cmq*Cm5l is sufficiently small. 
Notice that when [Cmg+Cmg,] + 2te/@ is greater than zero and if the non- 
linear term in the restoring moment can be ignored, equation (56) takes 
on the form 


'& - 286(1-ya*) + K = 0 (57) 


Equation (57) is then recognized as being the well-known Van Der Pol 
equation of nonlinear mechanics. It is evident that for small values 
of a the damping term is negative, leading to a divergent oscillation, 
whereas for larger a, the damping term is positive. A stable regime 
therefore will exist where a= l/u and oscillations of either large 
or small amplitude will converge to that regime. 


Thus, by the assumptions of an additional nonlinear pitching-moment 
variation caused by flow separation and a sufficiently small damping 
moment in the absence of separation, the observed phenomenon of a small- 
amplitude self-sustained oscillation can be hypothetically explained as 
being due to the destabilizing effect of the nonlinearity in the restor- 
ing moment on the damping moment. 


"It is recognized that usually the nonlinear effects of shock-wave 
boundary-layer interaction and flow separation are markedly reduced 
during nonsteady motions. However, for the low-frequency oscillations 
of the present tests (of the order of one cycle per 100 chord lengths 
of travel), it is believed that the effects of separation, though per- 
haps less severe than for the steady case, are nevertheless still 


present. 
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CONCLUSIONS 


The results of & theoretical and experimental investigation of the 
single-degree-of-freedom damping in pitch of a series of low-aspect-ratio 
wing-body combinations made at subsonic and supersonic speeds lead to the 
following conclusions: 


l. In the Mach number range 1.2 to 1.9, theoretical and experimen- 
tal values of the damping-in-pitch coefficient Cmq + Cm were in good 
agreement for triangular wings of aspect ratios 2, 3, and 4, and swept 
wing of aspect ratio 3. 


2. Theoretical predictions of the existence of ranges of center- 
of-gravity positions for a range of Mach numbers greater than 1.2 over 
which dynamic instability may be expected for the aspect ratio 4 trian- 
gular wing and the aspect ratio 3 swept and unswept wings were confirmed 
by the experimental results. 


3. The prediction by the theory of the beneficial effect on the 
damping in pitch at supersonic speeds of a reduction in aspect ratio was 
borne out by the results of experiments for the trianguler wings having 
aspect ratios 2, 3, and 4. In the Mach number range 1.2 to 1.9, the 
aspect ratio 2 wing had the largest damping moments throughout the range 
of Mach numbers and the aspect ratio l wing, the smallest. 


4. Experimental results for the stability boundaries in the super- 
sonic speed range of three wings of aspect ratio 3 having triangular, 
swept, and unswept plan forms confirmed the theoretical prediction 
regarding the relative magnitudes of the region of Mach number and 
center-of-gravity position in which dynamic instability could be experi- 
enced. In the Mach number range 1.2 to 1.9, no region of instability 
existed for the triangular wing. The region of instability for the 
unswept wing was considerably larger than for the swept wing. 


5. The effects of profile thickness on the damping in pitch at 
supersonic speeds of triangular wings of aspect ratio 4 and thickness 
ratios of 3 percent and 6 percent were found to be small and in the 
direction of decreasing stability with increasing thickness ratio. 


6. The occurrence at a Mach number of 0.9 of small-amplitude self- 
sustained oscillations of the triangular wing of aspect ratio 4 and the 
triangular and swept wings of aspect ratio 3 was attributed to a desta- 
bilizing effect on the damping moment of nonlinearities in the aero- 
dynamic restoring moment. 


Ames Aeronautical Laboratory 
National Advisory Committee for Aeronautics 


Moffett Field, Calif. 
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APPENDIX A 


TRANSFER OF AXES 


The problem to be solved may be posed ag: Given the damping 
coefficient Cmq + Cos and the restoring-moment coefficient Cm. about 


two axes of rotation, a known distance apart, find the damping coefficient 
about an arbitrary axis. 


For single-degree-of-freedom rotary oscillations about axis (1) the 
damping coefficient is written as 


= = = 2 
Cinq + ong, | = | omg T ong, T B3 | Cna | — Er. T "| — OX, SES (A1) 
1 O o 0 
where the axis (o) is chosen as & reference axis, | Cn + Ong]. ie known, 


and xX, is the known (nondimensionalized) distance of axis (1) ahead of 
axis (o). For pitching about axis (2), the damping coefficient is 


| na + 6 = | Ong + ong | + Ba Cna | = ža Cr, q * C | — 23e Op (A2) 


where | Ca 4 os and X2 are known. 


The relationship of the axes (1) and (2) to the reference axis (o) is 
shown in sketch (hh). Note that if the reference axis (o) is chosen to 


be coincident with the axis passing through +3 
the aerodynamic center, then the quantities * * 7 | 
t 
= = X 
Bac, and atal Om, | in equations (Al) K 
2 o (1) (2) (0) 


Sketch (hh) 
are known, the parameters Cig and Cm, may be determined from 
2 


and (A2) are zero. Since Cd and Cm, 
2 


C . 
"a *1 — * (A3) 
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With oT and om, | determined, it is now possible to solve equa- 


0 
tions (Al) and (A2) simultaneously. 


l 
w 


Let lone + Cm 


[——714 
Q 
B 
Q 
＋ 
e 
ao}. E 
H 
on 


2 = d 
Em 
a-b — e[X, — €,] + a[xi* — x47] 
Then E + 66 mild qc —————————— = e 
E " X, — 72 


(Al) 
| en = Ong, 3 2 — FI [-e] 4 dx = F 
o 


The damping coefficient about an arbitrary axis, with x referred to 
axis (c) and measured as positive ahead of axis (o), is therefore, 


| Cn + Ong, = f+ že- že- x (45) 
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Figure |. — Theoretical single - degrae -of - freedom short-period pitching stability boundaries 
for the two- dimensional wing at subsonic and supersonic speeds. 
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Figure 2.- Comparison of theoretical single -degree -of-freedom 
short-period pitching stability boundaries af supersonic 
speeds for three wings of aspect ratio 3, having 

friangular, swept, and rectangular plan forms. 
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Figure 3.— The effect of frequency on fhe theoretical single- 
degree — of - freedom short-period pitching stability 
boundaries at supersonic speeds for a triangular wing 

- of aspect ratio 4. 
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Figure 4.— The effect of frequency on the theoretical single- 
degree -of - freedom shor! -períod pitching stability 
boundaries at supersonic speeds for a rectangular 

wing of aspect ratio 4. 


NACA RM A5?2IOha 


Axis af "e 
: leading edge 


Reduced frequency, k 


Figure 5.— The variation with reduced frequency of fhe 
single- degree-of- freedom rotory damping — moment 
: coefficient for the two-dimensional wing at Mach 
numbers O, 0.50, and 0.80. Axis at leading edge. 


A 


LO Lf L2 L3 1.4 L5 J. L7 J. G S8 20 
Mach number, M, 


Figure 6.— The effect of thickness on the damping -íin-pifch coefficient Cm, for a two-dimensional 
wing having a symmetrical biconvex parabolic-arc section. Axis at leading edge. 
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Figure 7.— The effect of thickness on the damping -in-piftch coefficient Cm 7 + Cm 4 for a iwo- 
dimensional wing having a symmetrical biconvex parabolic—arc section. Axis at leading edge. 
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(c) A=4 triangular wing (3% biconvex circular - ore sections). 


Figure 8.— Continued. 
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Figure 8.— Conchided ; 
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Figure 9.— Experimental damping-in-pitch coefticients for the wing- body 
combination having a triangular wing of aspect ratio 2. 
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Figure 9.— Concluded. 
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Figure 10.— Experimental damping-in-pitch coefficients for the wing- 
body combination having a triangular wing of aspect ratio 3. 
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Figure 10.— Concluded. 
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Figure !l. — Experimental damping-in- pitch coefficients for the win g- 
body combination having a triangular wing of aspect ratio 4. 


I—— 


. NACA RM A552 TO 


Experiment : 
o R«133x10% 


4 4 
Axis t. 35€ 


Theory: —— 


Theory, wing + body + 
interference 


Nonlinear 
variation 
æ *5? fo /? 


6 8 LO L2 14 L6 L8 20 
Mach number, M, 


(b) Axis af O.35c. 


Figure ffl. —Confinued. 
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Figure !l.— Concluded. 
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Figure /2.— Experimental damping-in-piteh coefficients for the wing- 
body combination having a swept wing of aspect ratio 3. 
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Figure /2.— Concluded. 
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Figure /J.— Experimental damping-in-piteh coefficients for the wing- 
body combination having | an unswep! wing of aspect ratio 3. 
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Figure !3.— Continued . 
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Figure /3.— Concluded. 
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Figure 14.— Experimental pitching- moment -curve slopes for the wing- 
body combination having a friangular wing of aspect ratio 2. 
Axes of rotation af O.356 and O.45c. 
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Figure /5.— Experimental pitching - momen! -curve slopes for the wing- 
body combination having a triangular wing of aspect ratio 3. 
Axes of rotation at O. 35d and O.35c. 
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Figure 16.— Experimental pifchíng - moment! -curve slopes for the wing- 
body combination having a triangular wing of aspect ralio 4. 
Axes of rotation at O.256, 0.356, and O.45c. 
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Figure !7.— Experimental pitching moment -curve slopes for the wing- 
body combination having a swept wing of aspect ralio 3. 
Axes of rotation at | O.256 and O.35c. 
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Figure 18.— Experimental pifching -moment -curve slopes for the wing- 
body combination having on unswept wing of aspect ratio 3. 
Axes of rotation at 0.20¢, 0.35¢€, and 0.406, 
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Figure /9.— Experimental lift -curve slopes for the wing -body 
combination having a fríangular wing of aspect ratio 2. 
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Figure 20.—Experimental lift -curve slopes for the wing-body 
combination having a triangular wing of aspect ratio 3. 
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Figure 2l.—Experimental lift -curve slopes for the | wing-body 
combination having a triangular wing of aspect ratio 4. 
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Figure 22.—Experimental lift -curve slopes for the wing- body 
combination having a swept wing of aspect ratio 3. 
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Figure 23.— Experimental liff-curve slopes for the wing -body 
combination having an unswep! wing of aspect ratio 3. 
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Figure 24.— Comparison of experimental damping -ín-plIch coefficients for three wing -body 
combinations having friangular wings of aspect ratios 2, 3, and 4. Axes at 0.356. 
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Figure 25.— Comparison of theoretical and experimental single- 
degree -of- freedom short-period pitching stability 
boundaries for a wing-body combination having a 

triangular wing of aspect ratio 4. 
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Figure 26. — Comparison of magnitudes of experimental rotary damping- moments for fhree 
wing —body combinations having triangular wings of aspect ratios 2, 3, and 4. 
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Figure 27.— Comparison of experimental single- degree -of —- 
freedom Short period pitching stability boundaries for two 
wing — body combinations having swept and unswept wings 

of aspect ratio 3. 
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Figure 28.— Comparison of magnitudes of experimental rotary damping moments for fhree 
wing body combinations having triangular, swept, and unswept wings of aspect ratio 3. 
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Figure ES. Comparison of experimental damping-in-pitch coefficients for three wing - body 


combinations having friangular wings of aspect ratio 4 and thickness ratios of 3 and 6 percent. 
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Figure .30.~ Oscifograph record at Mach number 0.90 for the ung body combination 
having a triangular wing of aspect ratio 3. Axis at 0.25¢. 
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